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coefficients appearing in equations (5) 


Aj.Aj 

coefficients appearing in equations (94) 
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and (99), 

a 

length of panel 
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expressions defined by equations (44a) 
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expressions defined by equations (8) and 
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B 

set of bending modes and frequencies 
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width of panel 
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and appendix A) 


dj^j elements of beam frequency determinant obtained by approach 

of reference 7 (see eq. (102)) 


E£ Young's modulus of faces 

Ej,Fj analytic functions of 
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e 

distance to elastic axis o£ panel cross section measured 
from middle surface of core (see fig. 2 and appendix A) 

£ 

thickness of face sheets 
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Cy 

shear moduli of core in x- and y-directions, respectively 
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defined by equations (22c) 

K ) 

imaginary part of ( ) 

lo 

mass density moment of inertia of face sheets, per unit 
width, about elastic axis of panel cross section (see 
appendix A) 

i = 


i. j 

integers when used as subscripts 


partitioned matrices in equations (70) 

Nxb^ Nxa^ 

in-plane stress parameters, — — , — - — 

tt^Ds tt^Ds 

ky,ky 

Nyb^ Nya^ 

in-plane stress parameters, -4 — , -4 — 

ir^Ds tt^Ds 

Lij 

differential operators defined by equations (83) 

M 

Mach number 

m,n 

integers designating number of sinusoidal half-waves in x- 
and y-directions, respectively 

nTj = TTZj 


Nx’Ny 

intensity of in-plane force resultants (per unit width) 
parallel to x- and y-axes, respectively (positive in 
compression) 

*TJ 
1 — \ 
N 

II 

O 

characteristic equation 

P 

lateral loading on panel 

Qx>Qy 

intensity of internal shears (per unit width) acting in 
z-direction in cross sections originally parallel to yz- 
and xz-planes, respectively 

q 

dynamic pressure of airflow 

RC ) 

real part of ( ) 


X 



r,r ’ 

TT^Dg 

shear flexibility parameters, ^2j) ^ g^2Q 

s,s* 

defined by equations (22c) 

TS 

set of thickness -shear modes and frequencies 

TT 

set of thickness -twist modes and frequencies 

t 

time 

w 

transverse deflection of panel, measured in z-direction 

x,y,z 

orthogonal coordinates (see fig. 1) 


defined by equations (12) 

ttj = R(mj) 

roots of characteristic equations 

“x 

rotational contribution to panel slope in xz-plane (see 
appendixes A and B 

e = v4^ - 1 


« ^ 


S. = ICmj) 


Yx>Yy 

shear angles in xz- and yz-planes, respectively (see 

j • A ^ Qy 

appendix A) , 

Hx ^y 

C 

quantity defined by equation (18b) 

n 

panel length-width ratio, ^ 

A,X' 

j • * 2qb^ 2qa^ 

dynamic pressure parameters, ^ 

M 

/m \ 2 

Poisson's ratio of face sheets 


p = 


XI 



mass density of core 
mass density of face sheets 

mass per unit area of sandwich panel (see appendix A) 


Pc 

P£ 


Pm 


T = 


~s 


(see appendix A) 


$ = 





y 




x,x' 


X 


defined by equations (22c) 


• 4- 7T^ ^O 7T^ ^o 

rotary inertia parameters, , — ^ 

^ Pm ^ Pm 


function related to w by equation (91) 





co,oj' 

cr 

1,2 


- rxco^ . 

r- + n^r 

(zj/nn)^ 

panel circular frequency 

fundamental frequency of simply supported panel, rigid in 
shear, with an infinite length or width, respectively, 

7T^ 

/ Pm ’ ^ / Pm 
frequency parameters, ^ 

Subscripts 

critical (minimum) flutter value 
faces 1 and 2, respectively, when used with f 
properties in x- and y-directions, respectively 
Superscript 
complex conjugate 


Xll 



Notation 



3 ^ + ll. 

3x^ 9y^ Sy** 

A comma preceding a subscript denotes differentiation with respect 
the subscript. 




I . INTRODUCTION 


One of the basic structural elements used for exterior skin 
surfaces of flight vehicles is the flat rectangular panel. The flutter, 
or dynamic instability, of such panels when exposed to a supersonic 
flow has received considerable attention since the appearance of the 
pioneering works of Hedgepeth (ref. 1) and Movchan (ref. 2). 

A. Elastic and Inertia Considerations 
Most papers in the literature have dealt with homogeneous 
isotropic and orthotropic panels. However, high-efficiency structural 
configurations for supersonic vehicles may incorporate sandwich type 
panels having lightweight orthotropic cores such as honeycomb. Such 
panels often exhibit properties which are not accounted for by classi- 
cal plate theory. The primary difference between a homogeneous and a 
sandwich type panel is that the sandwich panel may experience signifi- 
cant transverse shear deformations due to the lightweight (low 
stiffness) cores that are often used. In homogeneous panels, this 
shearing deformation is usually negligible. 

The first analytic work dealing with the flutter of sandwich 
panels appeared in 1964 in a report by McElman (ref. 3). This analysis 
dealt with simply supported flat and cylindrically curved sandwich 
panels for which a two-term (mode) Galerkin solution was obtained. In 
the case of homogeneous panels, this type of solution was known to give 
qualitatively correct trends, but was not adequate (using only two 
modes) to predict reliable quantitative results (ref. 1) . One of the 
qualitative results of reference 3 was that the transverse shear 


flexibility of an isotropic core could significantly reduce the 
dynamic pressure required for flutter. 

In 1965 an exact closed form solution was published by Anderson 
(ref. 4) for the flutter of a flat rectangular, isotropic sandwich 
panel having simply supported or clamped leading and trailing edges and 
simply supported side edges. This investigation confirmed the signifi- 
cance of the isotropic core shear flexibility and showed that the two- 
mode Galerkin solution of reference 3 became progressively less accurate 
as the core shear flexibility increased. The theory developed in 
reference 4 was employed in reference 5 where numerical results for a 
wide range of panel parameters were presented for the simply supported 
edge condition. 

The numerical results obtained in references 4 and 5 indicated 
that the panel theory upon which they were based might be inadequate 
for certain ranges of parameters. These results showed that under 
certain conditions involving finite core shear stiffness, an increase 
in shear stiffness or boundary restraint would cause the dynamic pres- 
sure required for flutter to be reduced rather than increased. Because 
of these unexpected results, it was suggested that a more refined 
theory describing the panel motion might be necessary. For instance, 
one of the assumptions of the panel theory used in reference 4 is that 
the faces behave as membranes (i.e., the face sheet bending stiffness 
is neglected) . Consequently, the governing differential equations turn 
out to be sixth rather than eighth order, with a resulting loss of one 
boundary condition per panel edge. An additional refinement in the 
theory, suggested in reference 5, would be to account for the rotary 
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inertia. The inclusion of either the transverse shear flexibility or 
the rotary inertia in the analysis will lower the in-vacuo panel bend- 
ing frequencies. The simultaneous inclusion of both the shear flexi- 
bility and the rotary inertia, while still affecting the bending motion, 
will also admit two additional motions. These motions, which were 
described in reference 6 for a homogeneous isotropic plate, are termed 
the thickness-shear and thickness-twist modes. In general, both of 
these motions produce transverse deflections (perpendicular to the panel 
surface) and influence the aerodynamic loading. The question is whether 
these additional motions can significantly alter the flutter solution. 

The first consideration of the effects of the face sheet bending 
stiffness was presented by Smirnov in references 7 and 8 for a sandwich 
beam and panel, respectively. Unfortunately, the flutter analysis of 
the panel was based on an incomplete differential equation and the solu- 
tions presented are correct only for special cases. 

Flutter boundaries for sandwich panels having orthotropic shear and 
bending stiffnesses have been presented by Weidman (ref. 9) for a simply 
supported panel. As in reference 3, this analysis is based on a two- 
mode Galerkin solution and gives only qualitative trends. 

B. Aerodynamic a nd Damping Considerations 

As in references 1, 3, 4, and 5, the flutter analysis presented 
herein neglects all damping forces and is based on a simple two- 
dimensional static y inviscid approximation for the aerodynamic forces. 

In this approximation the aerodynamic loading is proportional to the 
instantaneous slope of the panel taken in the direction of the airflow. 
The analyses of references 2, 7, and 8 are based on a two-dimensional 
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quasi-steady aerodynamic approximation. This has the same mathematical 
form as the static approximation except for an additional term which is 
proportional to the panel's transverse velocity and represents an aero- 
dynamic damping force. Flutter boundaries based on these static and 
quasi- steady aerodynamic theories have recently been compared with 
results obtained from linearized three-dimensional, unsteady, inviscid 
(exact) aerodynamic theory (refs. 10 and 11). 

For unstressed homogeneous panels, all three theories predict 
nearly the same flutter dynamic pressure for Mach numbers greater than 
about 1.6 and for length-width ratios from 0 to at least 6. For Mach 
numbers less than 1.6 and for length-width ratios less than about 2, the 
"exact" theory predicts significantly lower flutter dynamic pressures 
than the other two theories. These lower values are due to a so-called 
single-degree-of -freedom instability and are very sensitive to struc- 
tural damping. (The other two theories do not predict this type of 
instability.) For these low Mach numbers all three theories tend to 
predict flutter boundaries that disagree with experimental results. 

This disagreement is apparently due, in part, to the stabilizing effect 
of the turbulent boundary- layer that occurs in this range of Mach num- 
bers and length-width ratios (refs. 12 and 13). For the higher Mach 
nimiber range (greater than 1.6) where the theoretical results agree, 
there is also fairly good agreement with experimental results (obtained 
at Mach numbers from 2 to 5 and length-width ratios from 0 to 10, 
ref. 11). For panels that are unstressed or are in tension, the results 
presented in references 10 and 11 also indicate that, for the higher 
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Mach number range, the flutter boundaries are insensitive to either 
structural or aerodynamic damping. 

For compressively stressed, unbuckled panels, the results of 
references 10 and 11 indicate that for Mach numbers greater than about 
1.6, the quasi- steady and exact theories predict flutter dynamic pres- 
sures that compare very satisfactorily . The static aerodynamics also 
gives essentially the same results except when compressive stresses 
cause the in-vacuo values of two bending frequencies to become equal . 
When this occurs, a zero-dynamic-pressure flutter point is predicted by 
the static aerodynamics. If the stress levels producing these equal 
frequencies are less than the buckling stress, then the quasi-steady and 
exact theories predict small but nonzero flutter values of dynamic 
pressure. (The zero points predicted by the static theory are prevented 
by the aerodynamic damping.) However, even these nonzero flutter 
dynamic pressures significantly underestimate experimental results 
(ref. 11). In fact, when an equal-frequency stress level coincides with 
the transition-point value of stress (buckling point of the flutter 
boundary) even the quasi-steady theory predicts zero-dynamic-pressure 
flutter (ref. 10, figs. 4 and 5). 

In reference 10 it is shown that reasonable quantitative agreement 
between theory and experiment can apparently be obtained near equal - 
frequency stress-level conditions, by the use of quasi-steady aerody- 
namics, if a structural damping term (g^) associated with the panel 
bending stress is included in the analysis. In this approach the zero- 
dynamic-pressure values for transition-point flutter do not occur. The 
reason for this behavior is not completely clear though since the 
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further addition of a structural damping term (g^^) associated with 
stresses produced by in-plane loads gives, for gg = flutter bound- 
aries which again have zero -dynamic -pres sure values for transition point 
flutter. 

To recapitulate, except for cases where compressive stresses cause 
two in-vacuo bending frequencies to have nearly the same value, the 
simple static aerodynamic approximation gives nearly the same results 
as the quasi- steady and ’^exact^^ aerodynamic theories for Mach numbers 
greater than about 1.6. It seems likely then that the results presented 
herein would not be significantly altered, in most cases, by the use of 
the more refined aerodynamics. For those cases of nearly equal in-vacuo 
bending frequencies, however, it would be appropriate to supplement the 
present analysis with a consideration of damping forces. 

C. Scope and Purpose of Investiga tion 

In this investigation an analysis is presented for the supersonic 
flutter behavior of flat rectangular, biaxially stressed sandwich 
panels. The mathematical model of the panel accounts for the bending 
stiffnesses of dissimilar face sheets, rotary inertia, and the ortho- 
tropic transverse shear stiffnesses (moduli) of cores such as honeycomb. 
In the absence of damping forces, for simply supported edges parallel to 
the airflow, and for the static aerodynamic loading approximation, an 
exact closed form solution to the problem is obtained for panels with 
either simply supported or clamped leading and trailing edges. The 
purpose of the investigation is to determine the effects and importance 
of the face bending stiffness, rotary inertia, and core shear stiffness 
orthotropy on the panel ^s flutter behavior. 
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II. ANALYSIS 


In this chapter the problem is defined and solved. The solution 
is found to take on different forms according to whether the panel has 
an orthotropic or isotropic core, or whether the length-width ratio is 
zero (beam behavior) . 

A. Panel Configurat ion and Differential Equations 
The panel configuration and coordinate system are shown in 
figure 1. The panel is flat and is of length a and width b. The 
panel is subjected to uniform, constant magnitude in-plane force resul- 
tants and N^. (positive in compression) . The in-plane shear force 
resultant, N^y, is assumed to be zero. At Mach number M, the flow is 
over one surface of the panel, is parallel to the x-axis, and is 
supersonic. Simply supported edges are assumed at y = 0 and y = b. 

The core is of uniform thickness and may possess orthotropic, 
transverse shear stiffnesses Dq ) . Its planes of orthotropic 

symmetry are parallel to the coordinate planes. The face sheets are 
isotropic but need not be identical (modulus, density, and thickness may 
differ but Poisson's ratio is assumed to be the same for both faces). 

The geometry of the core and face sheets is illustrated in figure 2. 

The quantities d^, d 2 , or o' locate the elastic axis. Stiffness and 
inertia properties are given in table I of appendix A. 

The extensional stiffness of the core is assumed to be much less 
than that of the face sheets (i.e., the core contribution to the overall 
panel bending stiffness is neglected) . Also, transverse shear strains 
in the faces and all normal strains in the z-direction are neglected. 
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The panel face sheets are considered to be thin plates rather than 
membranes; thus, their bending stiffness contributes to the equilibrium 
of forces in the z-direction (V^w term in eq. (la)). 

For the configuration described, an appropriate set of differential 
equations of motion can be obtained from either of the sandwich panel 
theories presented in references 14 and 15. This is done by incorpor- 
ating the effect of face sheet bending deformations in the sandwich 
panel theory of reference 14, or by incorporating the effect of ortho- 
tropic core shear moduli in the theory of reference 15. With these 
modifications, and with the addition of rotary inertia terms (see 
appendix A), both theories yield the following equations of motion. 

For equilibrium of forces in z-direction: 

”(Dfi Df2)V^w - “ ^y^,yy ^^xy^,xy '^Qx^'^Qy +p=0 

(la) 

For equilibrium of moments about y- and x-axes, respectively: 


-w - w 
,xyy ,xxx 




Qy 1 

W - W --^ + 77 

,yxx ,yyy Dg Dq^ 



(lb) 


(Ic) 


10 


I 



For equilibriiim of forces in x- and y-directions, respectively: 


-N^ 


+ N. 


xy 


= 0 


(Id) 


where 

Qx>Qy 

W 

P 


-Nv 


+ N. 


'^y,x 


= 0 


(le) 


bending stiffnesses of face sheets 1 and 2, respectively 
transverse shear forces, per unit width 
transverse deflection 
transverse pressure loading 

panel bending stiffness due to the extensional stiffnesses of 
the face sheets 


li Poisson's ratio of the faces (assumed equal for both faces) 

t time 

Iq mass density moment of inertia, per unit width, of the face 

sheets about the elastic axis 

In equations (Id) and (le) the in-plane inertia has been neglected. 

The last term in each of equations (lb) and (Ic) represents the rotary 
inertia moment about the elastic axis.^ It is assumed that the 
z-coordinates of the elastic axes and the centers of mass coincide. 

This condition is satisfied if the face sheets are of identical material, 
or if they have equal ratios of Young's modulus to density, Ef/p£. (If 
this condition is not satisfied, equations (lb) and (Ic) couple dynamically 
with equations (Id) and (le) , respectively (see appendix A).) 

recent flutter paper, in which the rotary inertia moments are 
incorrectly expressed, is discussed in appendix B. 
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The assumptions N^y = 0, Nx = constant, and Ny = constant satisfy 

the in-plane equilibrium equations, (Id) and (le), and eliminate the 

term 2N;^:yW from the lateral equilibrium equation (la) . The 
y ,xy 

unknowns, to be determined from equations (la, b, c) , are the lateral 
deflection w(x, y, t) and the two shear angles Qx(x,y^t)/DQ , 

Qy(x,y, t)/DQ^ (see eq. (Al)). The internal transverse shear force 
resultants Qx a-nd are depicted on the imagined panel cutout of 
figure 1. 

The transverse loading p(x, y, t) is comprised of the inertia 
force and the pressure due to supersonic flow. The aerodynamic loading 
is assumed to be given by two-dimensional static aerodynamics so that 
(ref. 1) 


P = 


-p w . . - w 
•^m ,tt g ,x 


( 2 ) 


where q is the dynamic pressure, g = - 1 , and is the panel’s 

mass per unit surface area. 


B. Boundary Conditions 

Equations (la, b, and c) constitute a system of eighth order 
partial differential equations in the spatial coordinates x and y. 
Corresponding to these equations are four boundary conditions which must 
be specified at each edge of the panel. Sets of boundary conditions 
consistent with the assumptions upon which equations (1) are based can 
be obtained from references 14 and 15. The boundary conditions used 
herein, in terms of the unknowns w, Qx/Pq ^ ^.nd Qy/Dq , are given below. 
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Simply Supported Edges 


w 


>x 


,xx 


Clamped Edges 


0 and a 

at X = 0 and a 


w = 0 

w = 0 

(3a) 

o 

If 

1 

-p- 

Qx ^ 

W - p; — = 0 

,x Dq^ 

(3b) 

w + yw =0 

,xx ,yy 

w =0 

(3c) 

W - = 0 

Qy n 

w - = 0 

(3d) 


For the simply supported case, equation (3b) specifies that no 
moment couple is produced by the face sheet axial loads (whose resultant 
is N^) . Equations (3c) specify that no net moment is acting on the 
individual face sheets. For the clamped condition, equations (3c) spec- 
ify a zero edge slope; hence, equations (3b) correspond to a zero shear 
angle in the x-direction. (Equations (3c) arise because the face 
sheets are treated as thin plates rather than membranes as was done in 
refs. 3, 4, 5, and 9.) 

Since equations (3a) are true at all points of the boundaries 
parallel to the y-axis, w ^ and ^ also zero at x = 0 and a. 
This, together with equations (3c), allows the boundary conditions to be 
written as follows. 
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Simply Supported Edges 
at X = 0 and a 


Clamped Edges 
at X = 0 and a 


w = 0 


w = 0 


(4a) 


Yx 





0 



(4b) 


w = 0 
,xx 


w 


= 0 


(4c) 


T 


y 



0 



(4d) 


Note that equations (4d) specify a zero shear angle in the 
y-direction. (See appendix A for the interpretation of the shear 
angles.) The boundary conditions at y = 0 and b are obtained from 

equations (4) by interchanging x and y. 

Along the edges x = 0 and a, the boundary conditions corresponding 
to equations (Id) and (le) require that (ref. 15): (1) either is a 

prescribed value or the translational displacement in the x-direction 
is zero; (2) either Nj^-y is a prescribed value or the translational 
displacement in the y-direction is zero. (The boundary conditions at 
y = 0 and b are obtained by interchanging x and y.) As noted earlier, 
the loads and Ny are assumed to be prescribed and the load Njjy is 

assumed to be zero. If and Ny are prescribed then the in-plane 

displacements will, in general, be nonzero. Uniform values of N^^^ and 
Ny could also be developed by uniform heating of the panel, in which 
case the in-place displacements could be zero. 
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C. _ Solution o^ Differential E quations for 
Simply Supported Streamwise Edges 
As in references 4 and 5 the solutions to equations (1) , which 
satisfy simple support boundary conditions at y = 0 and y = b, are 
sought in the form 


Z m (x/a) 
3 

^ in . (x/a) 


sin 


niry icot 
b ^ 


sin 


mry icot > 

b ^ i 


(5) 


Qy(x,y,t) 


E m. (x/a) 

C^e J 

3 


cos 


niry iwt | 

^ ' I 


The circular frequency is denoted by o), and n is an integer 
designating the number of half-sine waves that can form in the 
y-direction. 

The nTj are nonrepeated roots of the characteristic equation and 
are to be determined. As will be seen, the characteristic equation has 
eight roots, thus giving 24 values for the coefficients Aj , B ^ , and 
(j = 1 to 8) . However, the coefficients Aj , B j , and Cj are not all 
independent quantities since they are related by the three differential 
equations in w, Q^> and Q . Hence, one coefficient, say B • , can be 
expressed in terms of the other two (Aj and Cj) for each value of j. 
These eight coefficients (B^ to Bq) are then determined, within an 
arbitrary constant, by the eight boundary conditions (four at x = 0 
and four at x = a) . 

Substituting the above expressions for w, and into 
equations (1), and letting 
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C6) 


z, 

J IT 


results in the following set of three algebraic equations. 


— 



r 




an Zj 

-n 




0 

-Zj [n^(n2-Xur2)-Zj] n^-rx[zj - (i^)n^n^+n^Xar^j 

"j (^) 

< 

B, 

tt^Ds ^ 

>=< 

° f 

-n[n^(n2-Xw2)-z?] 



a^b 

ir^Ds ^ 


0 





L 

J 


L J 


C7) 


where 


an = ^ (n^n^ - zh^ + n^n^ky - k^z^ + z- 

j j ^ j 7T J 


(8) 


Thus, those values o£ zj which satisfy the determinant of the 3 by 3 
matrix will determine the nontrivial solutions to the differential 
equations . 

The dimensionless parameters appearing in equations (7) and (8) are 
defined as 


kx = 


= 


a 

b 

T 


Nxb^ 

tt^Ds 

ky 

_ ^y^ 

tt^D 

tt^Ds 


TT^D 

b^DQ^ 


~ b20( 

2qb3 

60s 

<^0 

TT^ 

■ b2, 

ir^ lo 

(J 

0) 

0)0 


Ds 

2 



(9) 
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The quantity n is the panel’s length-to-width ratio, and ky are 
measures of the in-plane loads, r^ and r^ reflect the core shear 
flexibilities, X characterizes the dynamic pressure, and X character- 
izes the rotary inertia. The ratio designated by x is a measure of 
the face sheet contribution to the transverse restoring forces and oT 
characterizes the panel frequency. 

The parameters appearing in equations (9) are all defined in terms 
of the panel width, b. For length-width ratios nil this is a conve- 
nient form, but for n < 1 it is more convenient to define all param- 
eters in terms of the panel length a. So, for n < 1 the following 
primed parameters will be used. 


k: = 


tt^Dc 


- ri^kx 


, Nya^ 


rv = 


tt^D 


a^D 


^ = 


Qx 


I tt^Ds 


""y a2o 


_ 


Qy 
2qa3 
BDs 


An' 


Wn = 



-I ^ 2 — 

i’ = — ^ = n 




Y t = li i = JL 


i2 P 


m 


an = -T 


[(nn )2 - Zj]^ + (nn)^ky - k^Zj + oj'^ - ^ Zj = n^a^ J 


(10) 


It is noted that n^Cn^ - Xco^) - z? appears as a common factor in 


the coefficients of Aj for both the second and third rows of 


17 


equation (7) . Also, n appears as a common factor in the coefficients 
of Cj for both the first and second rows of equation (7) . This 
enables equation (7) to be written as 





“ 


f \ 


' 

ah 

1 

1 

1 


1 -nn 
1 




0 

-z. ^z|-n2ri^+X'w'^+ J^yah^ 

1 

1 

1 

t 


1 

1 

1 

1 0 
1 
1 


tt3Ds 

► = ^ 

0 

0 

1 

1 

-nnxj 

1 

1 f 

1 7 • V • 


a3c- 


0 

) 

1 

( 

1 


TT^De 



(11) 


where 



The characteristic equation defining the quantities nTj = ttzj is 
obtained by equating the determinant of the square matrix in equa- 
tion (11) to zero. The form of the characteristic equation varies 
according to whether the core is ortho tropic, isotropic, or whether the 
length-width ratio is zero. 

1, Orthotropic Core (r^ ^ Vy) with r\ > 0. In this case the 
characteristic equation is 


an 




-[zj-(nn)^ + X'a)'^+(-^)rya^ J [z?yj - (nn) =0 


(13) 


t 4 ’ * 

Since a^^ contains z j , and and /j 
tion (13) is an eighth-order polynomial in 
nTj will appear in equations (5) . 


each contain 
Zj . Thus, 


2 

Zj, equa- 
eight values 


of 
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2. Isotropic Core 
expressions xj and yj 


(r^ = r ) with x\ > 0. For r^ = r' = r', 

V y 

are equal and equation (13) factors into 


the 


2„2 


[an - (z4 - n^n 


_,2 

Xo)' 


)(r'an + Zj - n^n^)]/! = 0 


(14a) 


The six roots of Zj given by the bracketed term depend on X' (dynamic 
pressure) . However, the two roots (to be designated zy and 23) given 
by 

y^ = 1 - r'[^(z? - n^n^) + X'aT'^] = 0 , j =7, 8 (14b) 


are independent of X * . 

3. Beam Behavior* (r\ = 0). Equation (13) also factors when n = 0 


(infinitely wide panel of 

length a) . 

However, 

it is 

to return to equation (11) 

and set n 

equal to 

zero t 


i "j 

1 

1 

1 

1 

1 

0 

. 2 . 1 2. 1+u » * ^ 

Zj(Zj+X'0) + 2 ^y^ll) 

1 » 2 l+]i 

' -X^+Z4 — ^ 

1 J J 2 

1 

(rx-ry) . 

0 

0 

! 0 
1 

1 

1 

1 



so that 


¥■* >=’”'1 


4 , , 2 _ ,2 X' 
-TZ-i -k^Z -+0)' - —r Zi 

J ^ J -n^ J 


-Zj (Zj+x'(jj'^) 


a^B^ 


tt^Ds 

a^Cj 


tt^D 


/ \ 
0 


>= < 


s/ 


0 

V / 


(15) 


-[l-rx(Zj+X'o)'^] 


= 0 (16) 


The determinant gives a sixth-order polynomial in Zj from which 
six roots (j = 1 to 6) can be determined. The bracketed term then gives 
Zj for j = 7 and 8, From equation (15) it is seen that Aj and Bj 
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are related but are independent of Cj . The 2x2 determinant 
corresponds to a "beam" solution since the same determinant is obtained 
by substituting the expressions 


w(x,t) 



Aje 


(x/a) 


iwt 

e 


Qj,(x,t) = 


6 


m.(x/a) . ^ 
„ 1 loot 

Bje e 


into the differential equations for the cylindrical bending of a beam 
(eqs. (la) and (lb) with all terms involving derivatives with respect to 
y set equal to zero) . 

D. In-vacuo Freque ncy Equations for All Edges Simply Supported 
For no airflow (A ' = 0) the differential equations and simple 
support boundary conditions are satisfied by 


w(x,y,t) = Ajnn 

sin 

mTTX 

a 

sin 

nTry 

b 

io)t 1 
e 



Qx(^>y>t) = Bmn 

cos 

mTTX 

a 

sin 

n-rry 

b 

icot 

e 

► 

(17) 

Qy(x,y,t) = Cmn 

sin 

m77X 

a 

cos 

nTry 

b 

iojt 

e 




where m and n are integers designating the number of half-sine waves 
that can form in the x- and y-directions, respectively. When equa- 
tions (17) are substituted into equations (1), with q = 0, the 
characteristic equations (13), (14), and (16) are again obtained except 
that A' = 0 and Zj = -m^. Thus, the characteristic equations become 
^frequency equations" from which can be obtained. 
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1. Orthotropio Core. The "frequency equation" form of 
equation (13) (orthotropic core) is cubic in co^. However, for X = 0 
(rotary inertia neglected) simplification occurs and the bending fre- 
quency of the m, n mode is given by 




= [(?) * (" * rir) - [(^) * "''‘y] 


where 






)rxr, 

[(m/T))2+n2]2 \ 

\ 2 

[rx(ra/n)2 +n2ry] 1 

1 

[-( 

¥)| 

5y(^) 



(18b) 


Equation (18a) includes the effects of length-width ratio, face sheet 
bending stiffness, orthotropic core shear flexibility, and in-plane 
stress. For t = 0 (face bending stiffness neglected) and rx = r^ 
(isotropic core) the above result reduces to equation (C4) of 
reference 5. 

2. Isotropio Core, For the isotropic core, equations (14a) and 

(14b) yield two frequency equations (if both the rotary inertia and 

shear flexibility terms X and r, respectively, are retained) . Replacing 
2 2 

Zj with -m in equation (14b) yields 


0)2 




(19) 


This equation gives the frequencies of the "thickness-twist" modes. 
These are free modes of vibration which, for a simply supported panel 
with an isotropic core, occur without any transverse deflection 
[w(x,y,t) = 0]. The frequencies (w, not oT) given by equation (19) are 
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independent of the in-plane loads (k^ and ky) , the face bending 
stiffness (t) and the mass per unit surface area (Pjj,) . The reason is 
that the transverse forces arising from the in-plane loads, face bend- 
ing deformations, and lateral inertia are zero when w(x,y, t) = 0 for 
all X and y (see eq. (la)). 

Except for differences in notation and interpretation of stiffnesses, 
equation (19) agrees with the second of equations (10) in reference 6. 
(The stiffnesses must be interpreted differently because reference 6 
deals with a homogeneous plate in a manner analogous to the Timoshenko 
beam theory.) Sketches of thickness-twist mode shapes for n = 1 and 
n = 0 are shown in figures 3(a) and 3(b), respectively. As pointed out 
in reference 6, the two shear angles Qx/Dq and Qy/Dq are so related as 
to cause the panel (n > 0) to twist about a normal to its surface. The 
equations for these mode shapes have been derived for A ^ 0 (eqs. (61) 
and (81)) and are found to be independent of the airflow. This is 
because w(x,y, t) = 0; hence, on the basis of the aerodynamic theory 
employed, no aerodynamic forces are produced. (This is not necessarily 
so for other boundary conditions or for a panel having an orthotropic 
core. ) 

The second in-vacuo frequency equation for the isotropic core is 
obtained from the bracketed portion of equation (14a) by setting A' =0 
and Zj = -m . This equation is quadratic in oT^ and yields 

(20a) 

where 


2[e^+(xp -^ (1+rg)] _ 

[l+(r+x) J+rX(Tj2_p-)] ± J Yl+(r-X)F-rX(TB^-fT)] 2+4x1' 
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(b) Thickness-twist modej section of infinitely 
wide panel, m = 2. 

Figure 3.- In- vacuo mode shapes for simply supported edges. 
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(20b) 



It is shown in appendix C that both sets of mode shapes corresponding 
to these two sets of frequencies involve bending and shearing deforma- 
tions for which w(x,y,t) 0. Thus, both sets of modes induce an 
aerodynamic loading that may lead to flutter. 

For convenience, the frequencies given by the smaller solution for 
(i.e., +y/ ) are referred to as the "bending" set of frequencies 

because for r = X = 0 (shear flexibility and rotary inertia neglected) 
they reduce to the frequencies of the pure bending solution given by 
classical plate theory. The larger solution for (i.e., -y/ ) 

gives the frequencies of the "thickness-shear" set of modes (described 
in ref . 6 for kx = ky = 0) . 

Additional information on these two sets of modes is given in 
appendix C where approximate bounds on the ratio of shear angle to total 
slope (Yx/'^,x) derived for a simply supported beam. The results 
indicate that the shear angle Yx larger than the slope w for 
the thickness- shear modes, whereas for the bending modes Yx less 
than w Y* The resulting forms for the bending and thickness-shear mode 
shapes are illustrated in figures 3(c) and 3(d), respectively. 

The in-vacuo frequency equations given by (18), (19), and (20) are 
extremely useful for the rapid estimation of frequency ranges of 
interest in the numerical computation of frequency loops (variation in 
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X with . These equations can also be used for the in-vacuo beam 
frequencies by rewriting them in terms of the primed parameters (see 
eqs. (10)) and letting n = 0. 

Equations (19) and (20) reveal that neither the thickness-twist nor 
the thickness-shear frequencies are predicted if either of the rotary 
inertia or shear flexibility terms (X and r, respectively) is neglected. 
The existence of the motions associated with these frequencies is due to 
the combined influence of the shear deformations and the inertia moments 
caused by rotational acceleration. 

E. Expansion of Characteris tic Equations 
in Powers of Zj 

In order to determine the roots zj from the characteristic 
equations it is necessary to expand these equations in powers of Zj . 
This is done in the following sections. 

1. Orthotvopio Core. This case is obtained by expanding 
equation (13) which results in the following eighth-order polynomial in 
Zj (j = 1 to 8) 

CgZ® + C0Zj + CsZj + C4Zj + C3Zj + C2Zj + ciz. + Cq = 0 (21) 

The coefficients Cj are given below in terms of the primed and 
unprimed parameters for n ^ 1 and 0 < n <. 1, respectively. 
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For nil 


C 8 = 
Ce = 

^5 = 

C 4 = 

C3 = 

C2 = 
Cl = 

Co = 

For 

C 8 = 
ce = 

C-5 = 
Clt = 

C3 = 
C2 == 

Cl = 
Co = 


-xr^ry /l-y\ 

^4 \ 2 ) 

[ts - r^ry(!-^)(kx - 2n2r) + (^)] 

(^)(^) 

(kx-2n2x) s-T (1+gr^) (1+hry) - (1+gry) + (<j)rxry-hry-n2rx) 


ns 


(22a) 


-n^|4>s+(kx-2n2x) (1+gr^) (1+hry) - |n2(l+gr^)+h(l+gry)+n2rxh^i^^|j 

-n^(l+grx) (1+hry) 

(1+grx) [(j) (1+hry) -n^h] 

0 < n <. 1 

-'Vy (^) 

(ki-2n2n2T)try(i^) 

{^) (^) 

(kx-2n2n2x) s.' -X (1+g 'r^) (1+h'ry) - (1+g 'ry) + c|) 'rxry-h'ry-n^n^r^) 

c./i: 


(t.'s' + (kx-2n2ri2x) (1+g’rx) (1+h'rp 

n2n2(i+gtj.^)+h' (l+g'ry) +n2n^rxh'^^-^^j| 

-( 1 +g'i’x) 

(l+g'r^^) [(j)* (1+h'ry) -n^n^h'] 


(22b) 
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The quantities appearing in 

the coefficients are given by 

(p = oT^+n^ky-n^T 

(Ji ' = to ' ^+n^n^ky-n‘^ri^T 


gi -X'to'^ 

h = n^ - Xu)^ 

h' = n^Ti^-X'to'^ 

s = (l+gry)rx+(l+hrx)ry^i^^ 

s' = (l+g'r^)r^+(l+h'rx)ry( j 


( 22 c) 


and the remaining parameters are defined by equations ( 9 ) and ( 10 ) . 

2 . Isotropic Core. For this case, six of the eight roots are 
obtained from the bracketed portion of equation ( 14 a) and must satisfy 
(j = 1 to 6 ) 


c 


SZj 


+ Cgz] + C^z^ + C^Z. + Cq 


= 0 


( 23 ) 


For n ^ 1 


C4 = 1 - rkjj + t(1 + 3 n^r - rXw^) 



C£ = n^{k3^-2n^(l+T)+r[n^(kx+ky)+^-3n‘*T]+X(J^(l-rk;x+2n^rT) } 
cj = n^(l + n^r - rxco^) 

Cq = -n^{ (w^ + n^ky - n‘+x) (1 + n^r - rXu^) - n^ (n^ - Xw^) } 


(24a) 
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For 


0 < n :! 1 


C 6 = -r'T 

C 4 = l-r’k^ + T(l + “ r 

C 3 = -r* 

^2 ” k^- 2 n^ri^ (l+T)+r * [n^ri^(kx+ky)+a7*^-3n^ri^T]+X*ar’^(l-r 'kx+ 2 n^r]^r*T) 

Cl = (1+n^ri^r* -r^X'co*^) 

cq =”{ (^' ^+n^n^ky-n^n^'r) (1+n^n^^ ' -r 'X 'Zo* -n^n^ (n^n^-X ’oT' } 

(24b) 

The other two roots are given by the term yj o£ equation (14b) 
and are 



zy = 
8 




Cl-r'x'oT’^) 


(25) 


3, Beam. This is a special case of either the orthotropic or 
isotropic panel solutions and is most directly obtained from equa- 
tion (16) • Six of the eight roots are determined from (j = 1 to 6) 


Cgz] + C 4 Zj + C 3 Z? + C 2 Z? + CjZj + Co = 0 (26) 


where 




^6 = 

C2 = kx + rxw'^ + X’ur'^Cl - r^kx)) 



C4 = 1 - r^k^ + t(1 - r^X’aT'^) 

Cl = (1 -’^xX'w'^)(^) 

► 

( 27 ) 

- = (^) 

Co = -rj^X'iir'^) 
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I III Hill nil 



The roots zy and zs are given by 



F . Relationships Between Aj , B j , and C j 

The coefficients Aj , Bj, and Cj appearing in equations (5) are 
not independent quantities. The relationships between them are deter- 
mined by equation (7) . As shown below, these relationships differ 
according to whether the panel has an orthotropic or isotropic core, or 
whether the length-width ratio (n) is zero. 

1. Panel. From the third of equation (11), which is just another 
form of equation (7) , the quantities B^ and Cj are related as 

Zj yjC^ - nnxjBj =0 (29) 

where Xj and yj are defined by equations (12) . For an isotropic core, 
Xj = yj so that equation (29) becomes 


[zjCj - nnBjlyj 


» t 


= 0 


(30) 


From equation (14b) it is seen that y] | t_ t = 0 only for j = 7 or 8 

j r^j-ry-r 

so that the bracketed term of equation (30) must be zero for values of 
j from 1 through 6. Thus, for the isotropic core. 



j = 1 6 


(31) 


where the Zj are the six roots of equation (23) . For the orthotropic 
core, equation (29) gives 
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Cj = f^4 B. j = 1 
Xj ^ 


where the zj are the eight roots of equation (21) . 

In terms of the primed parameters, the third of equation (7) is 


2 

-nn(-n^n^ + +z 


|)AJ . Bj 


Then, using equations (31) and (32), the relation between Aj and Bj is 


found to be, for the isotropic core 


Aj = r * + 


and for the orthotropic core 


l+r;(i 


n^n^ - Zj - X'w*^ 




)]f 'im-; 


j = 1 ^ 6 


(n^n^ - X'w'^ - Zj) 


©(; 


j = 1 8 


To determine Aj and in terms of B. when j = 7 or 8 (for the 
J J J 

isotropic core), the second of equation (11) is used 


(zj-n 


» [-xj » - 4)] Bj 


But for an isotropic core, r^ = r^ = r', and xj = = 0 for j = 7 

or 8. Since the coefficient of Aj is not zero in general, it follows 


Aj = 0 j = 7, 8 
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for the isotropic core. The first of equation (11) then gives 


C-i = B. 

J nn J 


j = 7, 8 


( 38 ) 


for the isotropic core. Except for the term X'oT' in Aj, the above 
relations between Aj, B j , and C ^ , for the isotropic core, are the same 
as in references 4 and 5. 

2. Beam. If equation (7) is rewritten in terms of the primed 
parameters (so that the column matrix does not contain the width b = ”) , 
the result for n = 0 is 


— 


1 n 





^iiln=o 


1 0 

1 


Aj 


'O' 

-Zj (Zj+X'Io'^) 

-[l-rx(Zj+X'u'^)] 

1 

1 0 

1 


a3Bj 


0 


TT^Dc 


0 

0 

1 v' 1 


II ^ 

a^Cj 



1 ^j^n=o 





(39) 


The third of equation (39) gives 


— ^ 


irSOe j‘n=0 


= 0 


But, from equation (16), 

only for j = 7 or 8 . Therefore, 


Cj-0 j=l ->6 (40) 

Thus, the only nonzero Cj are Cy and C 3 , and, as seen from equa- 
tion (39) , they are independent of Aj and Bj . 

It is only for values of j from 1 through 6 that the determinant 
of coefficients of Aj and Bj is equal to zero (giving eq. (26)). 
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( 41 ) 


Thus, Aj and Bj must be zero for j = 7 or 8. 


Aj = Bj = 0 


j = 7, 8 

For j = 1-6, the second of equation (39) gives 


■■ ■ (■; - 


j = 1 6 


(42) 


G. Satisfaction of Leading- and Trailing-Edge Boundary 

Conditions, and Corresponding Frequency Determinants 

Since the relationships between the coefficients Aj , Bj, and Cj 

have been determined, w, Q and Q can be expressed in terms of one set 

X y 

of coefficients. The boundary conditions can then be applied to obtain 
the frequency equations for the case of the orthotropic core, isotropic 
core, or beam. 

i. Orthotropic Core. When equations (32) and (35) are used, 
equations (5) become, for n ^ 1 and omitting cos|~^| 


m . (x/a) 


w 


fx) B. 

^ 7t3Ds bj Zj ^3 

^ j = l •’ ^ 


8 sr 


Qx(x) 


E m (x/a) 

e J Bj 

j = l 


m. (x/a) 


QyCx) = nn 


Zj ^1 

j=i 


(43) 


where x j , x j , yj, yj are given by equations (12) and a j , a^ , hy bj 
are defined as 


33 



(The primed quantities are to be used for n <. 1.) 

Application of the eight boundary conditions (eqs. (4)) to 
equations (43) gives the following set of eight homogeneous equations 


8 

^ ^ CijB-j =0 i = 1, 8 (45) 

where for n ^ 1 the coefficients c^.: are: 


Clamped Edges at Simply Supported Edges 

X = 0 and x = a at x = 0 and x = a 
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I£ 0 < ri < 1> a-j , bj, Xj , and are replaced with a j , 
respectively. If the determinant of the coefficients of 
equal to zero, the frequency equation for the panel with 
core is obtained. 


bj, x^. , and y j , 
Bj is set 
an orthotropic 


8x8 

lc..|=0 (47) 


2. Isotropic Core. When equations C^l) , (34), and (38) are used, 

sin f mry \ icjt 


equations (5) become, omitting the terms 


w ■ ;s5;2^(r * bj) 


w 


=E 




nij (x/a) 


Qx(^) = / ,e *' Bj 


j = l 


m . (x/a) 


Qy(x) = nn 


E m (x/a) 

B- + — / z^e ^ B- 

z • J nn / j J J 


J = l 


j=7 


(48) 


Application of the eight boundary conditions (eqs. (4)) to 
equations (48) gives the following set of eight homogeneous equations 


8 

c B. = 0 i = 1, 2, . . 8 (49) 

•1 — 1 J 


where the c^j are not the same as the c^j used in equation (45) . 
For nil, the coefficients are: 
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Simply Supported Edges at 
X = 0 and a 




3* Beam. In this case, w, Q^, and are independent of the 
y-coordinate, and equations (5) are replaced with 


w(x 


.t) = 


A^e 


io^t ^ 


ni , (x/a) . 

Qx(x>t) = ^ Bje ^ } (52) 

j 

(x/a) ^ 

Cje J 

j 

It is recalled from the discussion of equation (1^) that the above 
expressions for w and correspond to the *'beam’’ solution and lead to 
a sixth“Order polynomial in Zj = mj/Ti (given by eq. (26)). For 
w = Qx = 0 the above expression for when substituted into 

equations (1), yields the bracketed term of equation (16). The two 
roots of this term were called zy and zq (eq. (28)) and, as will be 
seen, correspond to the thickness-twist solution for an infinitely wide 
panel . 

When equations (40), (41), and (42) are used, the above expressions 
become (omitting e^^^) 

6 _ 

mj (x/a) 


w(x) = 


TT 


*■>3 


j=l 


QyCx) - i) 


J=7 


(53) 
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where 


b! = -(zf + X'w’ ) 

J J 

Application of the boundary conditions (eqs. (4)) gives two sets of 


homogeneous equations 
6 

SijBj = 0 


J = 1 


i 1^2^ • • 


(54a) 


O 


i = 7, 8 


(54b) 


J=7 


Of the eight roots for zj (given by eqs. (26) and (28)) the coeffici- 
ents Cj^j depend only on zj through Z5. These coefficients are given 
by: 



(55) 
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The coefficients Cj^j are the same for both sets of boundary conditions 
and are: 


C77 = C 7 e = 1 


(56) 


m- 


C07 = e 




^88 “ ® 


These coefficients depend only on zy and Z3. 

The frequency equation obtained by setting the determinant of the 
coefficients of Cj equal to zero is simply 


i-ij 


my 






I e ' e 

From equation (28) , iny and mQ are given by 


(57) 


m 


i = z. = ± 


1 - ryX'o)'^ 


j = 7, 8 


(58) 


For nonrepeated roots ^ equation (57) can be satisfied only if 


2 ^ 

1 - ryX'o)' < 0, so that 


mi 

Z) = -^ = ±i6, 

J TT 

Then equation (57) becomes 


/r^x'm' - 1 

<S = + / : > 0 


(^) ^ 


j = 7, 8 


-iTT6 imd T • X- n 

e - e = -2x sin iro = 0 


^For repeated roots (Z7 = zs = 0) the functional form of the 

0 fx/ a) 

solution changes to Qy(x) = [C7 + C8(x/a)]e ^ \ Application of 


Qy(0) = Qy(a) = 0 gives 
solution. 


1 0 
1 1 


^ 0, proving that Z7 = zs = 0 is not a 
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so that 


m = 1, 2 

j = 7, 8 


(59) 


m-j = ±i7r6 = ±im7r 


By squaring both sides of equation (58) , the use of equation (59) yields 

= -^(Jr + -- 2 ~ m = 1, 2, . . . (60) 

which is a special case (n = 0) of the thickness-twist frequency equa- 
tion. (See eq. (19).) Note that these frequencies are independent of 
the airflow. The mode shape corresponding to equation (60) is 


w(x) = Qx(x) = 0 


Qy (x) 
Dq 



(61) 


and is sketched in figure 3(b) for m = 2. 

The frequency equation for the bending and thickness-shear modes 
of the beam is obtained from equation (54a) and is 

6x6 

|£ij|=0 (62) 


The airflow does affect these frequencies since the zj, j = 1 ^ 6, do 
depend on X' . 

H. Uncoupled Solution for Simply 
Supported, Isotropic Panel 

The frequency determinants obtained for the panels (n > 0) determine 
the frequencies of the bending, thickness-shear, and thickness-twist 
modes. For the simply supported isotropic panel, the information 
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describing the thickness-twist motion uncouples from that describing the 
bending and thickness-shear motions. These uncoupled results are 
presented in the following sections. 

1. Frequenoy Determinant. The frequency determinant for the 
isotropic panel is given by equation (51) . For the simply supported 
edge condition this determinant can be simplified, by manipulating its 
rows and columns, to yield 
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(1 - rXu^)^ 

X 1 6 X 

6 1 X 

(e”8 - . 

3^7) 

= 0 


(63) 


(64) 


where |6 x 6 | is the 6 x 6 subdeterminant of equation (63). Note that 
it does not depend on the roots z-jim-j) and 23 (in's). 

From the first term of equation (63) it appears that the boundary 
conditions are satisfied for any values of the roots zj when the 

n 2 

frequency is such that 1 - rXu =l-r'x'w' =0. However, when 

n 

1 - rXu = 0 the characteristic equation (eq. (14a)), from which the 
roots Zj are determined, becomes 
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- n^rt^ + = 0 

which is seen to have two pairs of repeated roots (Z5 = Z7 = nn, 

Zg = Z3 = -nn) . The form of the solutions for w, and Q^. given by 

equations (5) is valid only for distinct roots of the characteristic 
equation. Thus, equation (64) is not valid at the frequency oT^ = 1/ (rx) . 
For < 1/ (rX) equation (64) is valid and becomes 

|6 X 6| =0 (65a) 

(e"*8 - e™^) = 0 (65b) 

The frequencies associated with equation (65b) are the thickness- 
twist frequencies. This is shown by considering equation (25) 


z 


j 



j = 7, 8 


( 66 ) 


For nonrepeated values of Zj and Zs, equation (65b) can be satisfied 
only if 


00 1 “ r'X'o)'^ r2 

n-^n^ + 




is a negative number (i.e., zy = -i 6, zq = +i6 ) . If such is the case, 
equation (65b) becomes 


+i7r6 -i7r6 

e - e 


2i sin 7 t 6 


0 


which is satisfied by 

6 — m • • • 


so 
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Zj = ±im 


j = 7, 8 
m = 1, 2, . 


Squaring both sides of equation (66) and using the above solution 
for Zj yields 

which is the thickness-twist frequency equation. The original deriva- 
tion of this equation (see eq. (19)) was based on X’ = 0 (no airflow), 
whereas the derivation above is valid for X* ^0. Note that these 
frequencies are independent of X (dynamic pressure) and therefore 
cannot coalesce with the bending or thickness-shear frequencies. (This 
is not necessarily so for other boundary conditions or for a panel 
having an orthotropic core.) 

Numerical results (ch. Ill) show that as X approaches zero the 
frequencies determined from equation (65a) approach the frequencies 
given by equation (20a) . Thus, equation (65a) gives the frequencies of 
the bending and thickness-shear motions. 

2. Mode Shapes. For all edges simply supported, equation (49) can 
be written as^ 


^For the first six columns of equation (68), the elements of 


column j are obtained from the corresponding elements of column one by 
replacing the subscript 1 with the subscript j . 
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j=l 2345 6 7 8 



Use of the relation between Aj and Bj for j = 1 6 (eq. (34)} 

enables the above equation to be expressed as 


( 68 ) 
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j = l 2345 6 7 8 
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r+(l/bi) • • 

r+(l/bg) 1 
1 

nn 
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nnke^^ 
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Zy 0 
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Zge 
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■ ■ r+(l/b6) 1 

nn 
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where k = Tr^Dg/nb^. Equation (69) is of the form 


6x6 6x2 6x1 


Kli j 

Ki2 


hi 


(o' 

1 

_K21 j 

K 22 

< 

,“jj 

^ , 

l4 


2x6 2x2 2x1 


so that 


[Kii]{Aj) + 

[Ki2]{Bj) = {0} 

(70a) 

[K2i]{Aj) + 

[K22]{Bj> = {0} 

(70b) 


Replacing b. with n^ - - (z^/h)^ (see eq. (44b)) and noting that 

J 

all elements of [K 12 ] are zero reduces equation (70a) to 
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The mode shape for w is then 


w(x,y,t) = 

mry 

sin 

E 

m . (x/a) 
Aje 3 


where, for o = 1 

(ref. 16 
1 

, p. 34) 
1 

3 = 1 

1 

1 

1 



e™3 

e^4 

e™3 

e”6 

Aj = i = (-l)^'^^c 
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^3 
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( 71 ) 


(72) 


(73) 


and so on, where c is an arbitrary constant. 

By setting the determinant of the square matrix in equation (71) to 
zero, one obtains, with manipulation of both rows and aotumns , 
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(74) 


Zge 


Except for the "missing” term (e™^ - e"^®) , equation (74) is the same 
frequency equation as given by equation (63) . 

Since equation (74) was obtained by altering columns (as well as 
rows) in the determinant of the coefficients of Aj , it is important to 
note that 
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(75) 


If one used equation (75), with an equal sign, the correct frequency 
equation for simple support boundary conditions would be obtained (see 
eq. (74)). However, the mode shape for w would be incorrect (see 
eq. (71)). 

The "missing" term (e™^ - e*"®) alluded to after equation (74) can 
be obtained from equation (70b), which is rewritten below. 
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(76) 


Consider the possibility that w = 0 for all values of the 
x-coordinate; then Aj = 0 for j = 1 -»■ 6 (in addition to Ay = A9 = 0 
by eq. (37)). Then, the above equation becomes simply 



giving a frequency equation 

ZyZs(e™® - e™^) = 0 


(77) 


(78) 


Although equation (78) is satisfied by zy or ze equal to zero, this 
solution cannot be considered valid since, by equation (25), zj = -zq 
which makes zy = Ze = 0 a repeated root. Thus, for distinct roots, equa- 
tion (78) reduces to equation (65b) which is satisfied by m^ g = limir, 
m = 1, 2, 3, . . . and again leads to the thickness -twist frequency 
equation (eq. (67)). 

Because Aj = 0 for j = 1 -> 6, equations (34) and (31) show that 
Bj = Cj = 0 for j = 1 6. The mode shapes given by equations (48) then 

reduce to 
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w(x) = 0 

Qx(x) = t 

Q„(x) . t is. 

nn ^ nn ° 

Since zy = -Z3 ^ 0, equation (77) gives 

By - Be = 0 (80) 

The use of equations (5), (79), and (80) together with in = ±im7r, 

7,8 

m = 1, 2, 3, . . •, then yields 
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(81) 


as the mode shape for the thickness-twist motion. (See fig. 3(a).) 

I. Single Differential Equation Approaches 
and Comparison with Smirnovas Results 
It is seen from equations (21), (23), and (25) that the 
characteristic equations for both the orthotropic and isotropic cores 
correspond to an eighth-order system of differential equations. In 
reference 8, Smirnov obtained flutter solutions for sandwich panels 
having isotropio cores; these solutions were derived from a sixth-, 
rather than an eighth-order differential equation. The sixth-order 
equation is an incomplete description of the panel motion, but for 
simply supported edges it does happen to lead to the correct bending and 
thickness-shear solutions. This occurrence is related to the uncoupled 
frequency determinant given by equation (63) . The connection between 
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the sixth-order solution for the simply supported, isotropic panel and 
the uncoupled eighth-order solution given by equation (63) is brought 
out in the following sections. 

1. Uncoupled Solution, In general, the solutions for w, Q^, and 
Qy. depend on all eight roots (irT) of the characteristic equation (see 
discussion after eq. (5)). For a panel having an isotropic core, how- 
ever, equation (37) shows that Kj = Kq = 0 even though Bj and Cj are 
all nonzero. Thus, in this special case, the solution for w depends 
on only six of the eight roots while the complete solution (i.e., for 
w, and Qy) still depends on all eight roots (as shown by eqs. (48)). 

With A 7 = As = 0, half the elements in columns 7 and 8 of the 
frequency determinant are zero (Fjj = ^ 2 ^ = ^ 5 ^ = = 0 in eqs. (50), 

j = 7 and 8) . Furthermore, for simply supported edges the elements 'Cgj 
and of columns 7 and 8 are proportional to elements and '^ 3 j , 

respectively. This leads to the uncoupled frequency determinant 
(eq. (63)) wherein m 7 and in’s are associated with the thickness-twist 
frequencies and m 2 through mg are associated with the bending and 
thickness- shear frequencies. 

The uncoupling of the thickness-twist solution from the bending and 
thickness-shear solution is seen to be caused by two conditions. The 
condition A 7 = Ag = 0 is a result of the factorable characteristic 
equation, (14a), and is due to the isotropic core. The proportionality 
of elements in the frequency determinant is due to the simple support 
boundary conditions . A change in either of these conditions prevents 
the solution from becoming uncoupled. Thus, a clamped isotropic panel 
will not have uncoupled solutions, nor will a simply supported 
orthotropic panel . 
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The fact that A7 = Ab = 0 arises from the factorable character- 
istic equation suggests that the eighth-order system of differential 
equations given by (la), (lb), and (Ic) also has a factorable form. 
This is shown next. 

2. Eighth-Order Differential Equation. The three differential 
equations (la), (lb), and (Ic) can be written in the form 
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and p is given by equation (2) . 


(82) 


(83) 
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is obtained by 


A single differential equation in w, or 

expanding the determinant of the matrix treating the operational 

coefficients as algebraic quantities (ref. 16, p. 117). The general 

result is an eighth-order differential operator, D 3 , operating on w, 

Qjj, or Q . For an isotropic core (Dq = Dq = Dq) the operators L 23 
y X y ' 

and L 32 are equal and the operator Ds factors into the product 

D 3 = D 2 D 0 (84) 

where D 2 and D 0 are the second- and sixth-order differential operators 



(85a) 

(85b) 


When the operators D 2 and D 0 are applied to any one of 

equations (5) and the results set equal to zero, equations (25) and (23) 

are reobtained, respectively, giving the e-lght roots z j . This approach 

has no particular advantage over the method of working directly with 

equations (82) since the three simultaneous equations must still be used 

to determine the relations between A.:, Bj, and C.. The formal 

J J’ j 

expansion of the determinant of [Lj^j], although lengthy,** does give the 
complete ei-ghth- order differential equation. If, however, the two shear 
angles (Qx/Dq^j Qy/Dq^.) eliminated from equations (82) by 


**This expansion is given in reference 14 for the case of an 
orthotropic panel, but with the face sheet bending stiffness and 
rotary inertia neglected. 


52 



differentiation, only a sixth-order equation in w is obtained. This 
is shown in the next section. 

3. Sixth-Order Differential Equation. For an isotropic core 
(Dq^ = Dq = Dq) , the sixth-order differential equation in w is 
obtained by first differentiating equation (lb) with respect to x, 
differentiating equation (Ic) with respect to y, and adding the 
resulting equations together to obtain 

'■o 3^ 


where 




( 86 ) 


(87) 


Equation (la) can be written as 


(Df^ + Df^)v'+w + Nxw^xx - p = 0 


( 88 ) 


by noting the identity 


Dn($ - v2w) = -CQx^x J 


>y 


The solution for ^ in terms of w is easily obtained from equa- 
tion (88) and when substituted into equation (86) leads to a sixth-order 
differential equation in w 

‘>s[i - 

- lo ^ v2w - ^ [(Df ^ + Df 2 ) v^w + Nx ^ v2w + Ny ^ v2„ _ v2pl = o (89) 
3t^ Q L -J 
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Equation (89) corresponds to DgCw) = 0 (eq. (85b)) since when p 
and w are replaced with equations (2) and (5) , respectively, equa- 
tion (23), giving six of the roots z j , is again obtained. The factor 
D 2 (eq. (85a)) is missing from equation (89) and hence the roots zy 
and zq, given by equation (25), are not obtained. Also, since equa- 
tion (89) is a sixth-order equation, only three boundary conditions per 
edge can be specified. Thus, the sixth-order differential equation is 
seen to be incomplete, and solutions derived from it alone are, in 
general , incorrect . 

4. Solutions Based on Sixth-Order Equation. For simply supported 
edges the incomplete sixth-order system does happen to lead to the 
correct frequency equation for the bending and thickness-shear modes. 
This can be shown by introducing a function X defined by^ 

$ = v^x (90) 


(The function X is used to express the boundary conditions in terms 
of a single variable.) Comparison of this result with equation (86) 
shows that 


w 


lo 32 
^ Dq at2 



If equation (2) is written as 


(91) 


S2w _ 

P = -Pm ^ - q 


( 92 ) 


^The symbol X corresponds to Smirnov's X (ref. 8) and is 
distinct from the X used herein for the rotary inertia parameter. 
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where q is the aerodynamic loading, equations (88), (90), (91), and 
(92) can be combined to give a sixth-order differential equation in 


1 + T - T — 
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v‘*x + (Nx -^ + Ny 
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When the rotary inertia terms [Io(^^/^t^)] are dropped, and x is 
neglected with respect to one, equation (93) reduces to the differential 
equation used by Smirnov (ref. 8). Replacing q with the aerodynamic 
loading term [ (-2q/g) (9w/9x) ] and setting 


X 
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— (x/a) . nuy io)t 

A-:e sin —r^ e 

1 b 

j = l 



(94) 


in equation (93) , once more leads to equation (23) for six of the roots 
z j ; again, zy and zq are not obtained. 

If, for simply supported edges, boundary conditions (4b) and (4d) 
are combined (so as to have three rather than four boundary conditions) , 
the following modified boundary equations are obtained. 


at X = 0 and a 


at y = 0 and b 


w = 0 


w = 0 


'^,xx = 0 


W^yy 0 


— CQx.X Q,. J = 0 ^ (Qy^y + = 0 


Dq '•^X,X Hy^y 


(95a) 

(95b) 

(95c) 


55 



These thvee boundary conditions per edge are satisfied along the panel 
boundary provided the function X is required to satisfy 

X = V^jT = v'+X = 0 (96) 

along the boundary . ® 


Applying equations (96) to equation (94) yields the following 
relationships between the six Aj . 
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Note that the determinant of the square matrix in equation (97) is 
identical to the uncoupled 6x6 subdeterminant of equation (63) . 

The relation between Aj (corresponding to )T) and Aj 
(corresponding to w) is determined by equation (91) and is found to be 

Aj = (1 + rbj)Aj j = 1 ^ 6 (98) 


^Equations (95c) are satisfied because of the identity 



where equation (44b) has been used. By employing equation (98), it is 
possible to manipulate the rows of equation (97) so as to reobtain 
equation (71) which was originally derived from the complete eighth- 
order system. 

Equation (97) depends on only six roots, and these are the same 
roots as are given by equation (23). Thus, the sixth-order equation, 
used with the modified simple support boundary conditions, does predict 
the correct solution for the bending and thickness-shear motion. (Note 
that the sixth-order equation does not predict the thickness-twist 
solution.) For other boundary conditions, however, (several of which 
are considered in ref. 8) the bending and thickness-shear motions will 
not, in general, be independent of the other two roots given by the 
second-order equation (i.e., eq. (85a)). In such cases the sixth-order 
system will not give the correct frequency equation. In reference 17, 
several other instances are noted where the assumption of core isotropy 
has led to the erroneous use of differential equations that are two 
orders less than is appropriate. 

Another instance in which the characteristic equation factors is 
when n = 0 (beam behavior) . In this case Aj and Bj are related but 
are independent of Cj (eqs. (40), (41), and (42)). This reflects the 
fact that the ’^cylindrical bending” of the sandwich beam is governed by 
a sixth-order rather than an eighth-order differential equation (see 
discussion after eq. (16)). 

The bending and thickness-shear frequency equation for the simply 
supported beam (n = 0) is given by equations (55) and (62) as 
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By manipulation of the rows, this determinant can be reduced to the 
determinant obtained from equation (97). Thus, the sixth-order system 
correctly predicts the bending and thickness-shear frequencies for both 
the simply supported beam and the panel. 

For the clamped beam (n = 0) the sixth-order differential equation 
also predicts the correct frequency equation. In this case 

X = i) C 99 - 

j=l 


and application of 


/i + £ 2 . AL\ 

\ Dq 3t;2 Dq dx^ / 


dx 


d^X 

dx^ 


= 0 


yields, 2 


( 100 ) 




i = 1, 2, . . . , 6 


1 = 1 


(101) 


where 


d. . = 1 - r 'X 'o) ' -r'z^ 

11 1 


d„ . = z . 
3l 1 


o 9 ni- ^ 

d^j = (1 - r 'X'cu' - r'Zj)e J 


m- 

d^. =z.e 1 


d^. = z? 
5l 1 


Q m j; 

d, . = z?e 3 
^1 1 


( 102 ) 


"^Except for the term (Iq/Dq) (3^/9t^) (rotary inertia) these 


boundary conditions are the same as used by Smirnov in reference 7. 



Hence, the corresponding frequency equation predicted by the sixth-order 
system is 

6x6 

|d.j| = 0 (103;) 

From the eighth-order system the bending and thickness-shear 
frequency determinant for the beam solution is 

6x6 

le.^l = 0 (104) 

(See eqs . (54a), (55), and (62).) The elements Cj^j are given by 
equations (55) and depend on only six of the eight roots for . For 
the clamped beam these elements are seen to be different from the 
elements dj^j predicted by the sixth-order system (eqs. (102)). By 
using equation (98) and by manipulating the rows of equation (104), 
however, it is possible to reduce the elements of 1 ^ | to the corre- 
sponding elements of [dj^j | . Thus, for the beam (n = 0) , the sixth-order 
system predicts the same clamped bending and thickness-shear frequency 
solution as the eighth-order system. (Note that the sixth-order system 
does not give the thickness-twist solution.) For n > 0 (panel rather 
than beam) the clamped eighth-order solution depends on all eight roots 
of the corresponding characteristic equation (see eqs. (49), (50), and 
(51)); in this case the sixth-order system does not adequately describe 
the motion. 

The above results show that for the simply supported isotropic 
panel and the simply supported or clamped beam, the solution from the 
eighth-order system of differential equations uncouples and can be 
expressed as a sixth- plus a second-order system. The solution is 
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uncoupled in the sense that the sixth-order system predicts the bending 
and thickness -shear frequencies while the second-order system indepen- 
dently predicts the thickness-twist frequencies. In the approach used 
by Smirnov a single differential equation in w corresponding to the 
sixth-order system only is obtained. This approach is therefore not 
applicable to the general case where the eighth-order system does not 
uncouple. 
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III. COMPUTATIONAL PROCEDURE FOR 
OBTAINING NUMER ICAL RESULTS 

For the type of aerodynamic-force approximation used herein, pairs 
of panel frequencies coalesce at certain dynamic pressures, thus form- 
ing ^frequency loops” in the X, o) plane. (See fig. 4.) Below the 
peak of a particular loop the two frequencies forming the loop are real 
and distinct. In this case the motion of the corresponding modes is 
stable. At the peak (i.e., where 9X/3o) = 0) the two frequencies have 
coalesced; they are still real at this point but are on the verge of 
becoming complex. Above the peak the pair of frequencies exist as com- 
plex quantities and produce an oscillatory diverging motion termed 
flutter (ref. 18, p. 245). The value of X corresponding to the peak 
of the lowest loop corresponds to the critical (minimum) flutter 
condition and is designated X^^. (See fig. 4.) 

Points on the frequency loops are combinations of X and 15^ for 
which the roots Zj of the characteristic equation also cause the 
corresponding frequency determinant to equal zero. A practical proce- 

O 

dure for determining the frequency loops is to fix oj and vary X 

between prescribed limits, computing first the roots Zj and then the 

frequency determinant D for each A. These steps can be repeated 

until a value of A is found that differs from the correct value (i.e., 

D = 0) by only an acceptable amount. (An allowable error in A of 

0.01 percent was accepted in calculating the numerical results presented 

herein.) By incrementing and repeating the above process one can 

systematically obtain the frequency loops. Of course, the same proce- 

2 

dure can be used by holding A fixed and iterating to solutions for m . 
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A. Roots of the Characteristic Equations 

In references 4 and 5 the characteristic equation was a 
fourth-order polynomial whose roots could be computed analytically. The 
characteristic equations obtained herein are eighth- and sixth-order 
polynomials (eqs. (21), (23), and (26)) and their roots were computed 
numerically, using Bairstow’s method as programmed in reference 19. As 
a check on the accuracy of the computations, the coefficients of the 
polynomial as computed from the roots were compared with the input 
coefficients; the coefficients generally agreed to 12 digits or better. 
(All computations were performed in double precision arithmetic.) 

There is no guarantee that the array of roots returned from the 
root finding subroutine will be ordered in a consistent, systematic 
manner. For example, when X = Xi the two roots stored in z(l) and 
z(2) may be the most negative and the second most negative roots, 
respectively. On the next trial value of A, X = + AX, the roots 

stored in z(l) and z(2) may be the second most negative and the most 
negative roots, respectively, with the remaining roots ordered in the 
same sequence as for A = A^ (see fig. 5). Computing the frequency 
determinants D(A) and D(A + AA) with the names (i.e., and Z£) of 
the two most negative roots interchanged when A = A2 + AA is equiva- 
lent to interchanging columns 1 and 2 of D(Ai + AA), hence producing a 
false change of sign. Since a sign change indicates the existence of a 
solution in the interval AA, a consistent, systematic procedure for 
ordering the roots must be incorporated in the programming to prevent 
the occurrence of false changes in the sign of the frequency determinant. 
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Figure 5.- Schematic of characteristic equation for two 
values of X with fixed 
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The approach used herein was to identify the real and complex roots and 
then to sort and name them according to the magnitude and sign of their 
real and imaginary parts. 

B. Frequen cy Determinant 

Some of the roots zj will usually be complex, thus causing 
elements of columns j of the frequency determinant to be complex also. 
Since the coefficients appearing in the characteristic equation 
P(zj) = 0 are real, however, complex roots will occur in complex con- 
jugate pairs. This allows the frequency determinant to be computed from 
a matrix of veaZ numbers. To show this, consider the 4x4 determinant 


El 

E2 

Es 

E4 


E2e”'2 

E3e“3 

E 46 

El 

F2 

F3 

F4 


F 26™2 

F3e“3 

Fi^e 


(105) 


where Ej = E(zj) and Fj = F(zj) are analytic functions of 
Zj = ~ (otj/ 7 T) + i(6j/TT) that are real when Zj is real. 

(Eq. (105) is of the same form as the frequency determinants derived 
herein.) Then if z\ and Z 2 are complex conjugates, let 


7 TZi - ai + i^i = m^ 


TTZ^ = 'n‘Z2 = a-j^ - i6^ = m^ 


so D becomes 
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I 



E(zi) 

E(z*) 

E 3 

E 4 

E(zi)e“l 

E(z*)e®l 

E3e®^3 

E 4 e”“* 

F(Zj) 

F(zp 

F 3 

F 4 

F(zi)e“l 

F(z*)e®^ 

F3e"3 

F 4 e®’‘* 


By the principle of reflection (ref. 20), E(z*) = [E(zx)]*, and so on, 
giving 


D 


E(zi) [E(zi)]* Eg E4 

E(zi)e®l [E(zi)e”'l]* Eae™^ Eite®** 

F(zi) [F(zi)]* F 3 F 4 

F(zi)e"*l [F(zi)e"'l]* Fae’"^ F4e"“* 


If column 2 is added to column 1, the quantity 2 factored from the 
new colimin 1 , and the resulting column 1 subtracted from column 2 , the 
above expression for D simplifies to 


R[Ei] I[Ei] E3 E4 

R[Eie^^l]e“^ I [Eie^'^l] e“l Eae®^ E4e™'+ 

D = -2i (106) 

R[Fx] I[Fi] F 3 F 4 

R[Fie^^l]e“l I[Fie^^l]e“l Fae"*^ F4e"'‘* 

where R and I denote the real and imaginary parts, respectively, of 
thfe bracketed expressions and where mi = + i 6 x has been used. 

Note that D is purely vmaginavy when 23 and 24 are real. If 23 and 
Z 4 are also a complex conjugate pair (7723 = m 3 = U3 + i63, 24 = 2p, 
then columns 3 and 4 can be manipulated in the same manner to give 
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R[Ei] 

I[El] 

R[E3] 

I[E3] 



R[E36^'^3]e“3 

I[E3e^'^3]e' 

R(Fl] 

I[Fi] 

R[F3] 

I[F3] 



R[F36^^3]e“3 

I[F3e^^3]e' 


which is a veal number. 

The above procedure applies to an nxn determinant whose elements 
are of the form of equation (105). Thus, the frequency determinants 
considered herein can be computed from determinants whose elements are 
real rather than complex. This result is useful since it gives one 
access to a wider range of available methods for evaluating the 
determinants . 

a . m4 

The terms e ^ and e ^ appearing in the frequency determinants 
can cause a computer overflow condition for large positive exponents 
(i.e., e^j may be too large a number for the computer to represent). 

A simple way to avoid this problem is to write the elements of the 
determinants so that no positive real exponents occur. For example, 
consider equation (106) with ai > 0, Tn 3 < 0, and m 4 - > 0. Then, 

equation (106) can be written 



R[Ei] 

I[Ei] 

E3 



R[Eie^'^M 

I[Eie^^l] 

E3e“3-«l 

E 4 

2 i 

R[Fi] 

I[Fi] 
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R[Fie^^l] 
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F3e“3-“i 

F 4 


D = A 
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where A is the determinant of the 4x4 matrix in the above equation. 
Since it is only the sign of the frequency determinant that is of inter- 
est (as D passes through zero), it is only necessary to compute A, 
which contains no elements having e raised to a positive real power. 

It should be noted that some of the elements of the frequency 
determinants derived herein involve terms where bj , y j , and appear 
in the denominator. If one of these quantities passes through zero (as 
X or ( 0 ^ is being incremented) the frequency determinant D will change 
sign by passing through plus and minus infinity, rather than by passing 
through zero. Thus, the form of the frequency determinant should be 
modified so as to prevent D from becoming infinite. This can be done 
by multiplying the elements of column j by the appropriate products 
of b j , y j , and Zj so that none of these three quantities appear in the 
denominator of any element. 
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IV. xNUME RICAL RESULTS AND DISCUSSION 
The procedure discussed in chapter III has been used to compute 
frequency loops and flutter boundaries for a selected range of param- 
eters. These results, which are discussed in the following sections, 
indicate the effects and importance of the face sheet bending stiff- 
ness, rotary inertia, and orthotropic core shear stiffnesses. 

A. Effects of F ace Sheet Ben di ng Stiffness 
The influence of the two face sheet bending stiffnesses Df^ and 
Df^ is felt through the parameter t = (Df^ + Df^)/Ds. If both face 
sheets are identical this parameter becomes t = (1/3) { (f/c) / [1+ (f/c) ] }^ 
which depends only on the face sheet to core thickness ratio (f = face 
thickness and c = core thickness) . 

In reference 4, flutter boundaries were calculated for sandwich 
panels having clamped or simply supported leading and trailing edges and 
simply supported side edges. These results, which were from an analysis 
that neglected the face sheet bending stiffness, predicted that for both 
a square panel (n = 1) and an infinitely wide panel (n = 0, i.e., a 
beam) the clamped edge condition gives a lower flutter speed than the 
simply supported edge condition when the shear flexibility parameter, 
r*, exceeds a value of about 0.2.^ 

®For an isotropic core and identical face sheets 



where Ef and are Young’s and the shear modulus of face and core, 
respectively. 
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That an increase in boundary restraint could lower the flutter 
speed was an unexpected result. Indeed, Smirnov's flutter analysis of 
the sandwich beam (ref. 7) showed that for "large" values of f/c and r' 
(f/c = 0.632, r' = TT^) the clamped edge condition does give a higher 
flutter speed than the simply supported edge condition. Since Smirnov 
did not present results corresponding to smaller values of f/c, it was 
not possible to determine whether the difference in the effect of the 
boundary conditions predicted by reference 7 (where t > 0) and refer- 
ence 4 (where t = 0) was due to the "largeness" of the f/c chosen by 
Smirnov (ref. 7), or whether it was the neglect of x in reference 4. 
Neglecting x lowers the order of the beam differential equations from 
sixth to fourth order (see eq. (26)) and raises the possibility that 

lim A^^(f/c) ^ ^cr^°^ (refs. 21 and 22).^ 
f/c-s-0 

To investigate this possibility, flutter boundaries were computed 
from the solutions to the sixth-order beam equations derived herein for 
clamped and simply supported edges, and for f/c ranging from 0.01 to 

0. 632. These results are shown in figure 6 where the variation in 

with f/c is presented for r* = 0, 0.4, and 2.0 (with = 0^ 

1. e., no in-plane load, and rotary inertia neglected). The curves 
labeled r’ = 0 correspond to panels having an infinite shear modulus, 
and in this case the only effect of x is to increase the values of 
A^r computed from the x = 0 analysis by the factor (1 + x) (see 

eq. (26)). (Note that Ds(l + x) equals the total panel bending stiff- 

ness.) For r’ = 0.4, the clamped beam has a lower flutter boundary than 
^By (f/c) ^ 0 it is understood that (f/c) goes to zero only in the 
expression for x. 
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the simply supported beam until f/c 0.53 where the boundaries cross. 

The same situation exists for r' = 2.0 except that the boundaries cross 

at a smaller value of f/c. In all cases lim A^j.(f/c) = Aqj.(0) 

f/c->-0 

indicating that the flutter solutions to the fourth-order beam equation 
presented in reference 4 correctly describe the beam bending behavior 
for sufficiently small values of f/c. For values of f/c greater than 
zero the face bending stiffness is seen to increase in importance as the 
shear flexibility parameter (r') increases. 

The behavior of the frequency loops for four points on the r' =2.0 
flutter boundaries of figure 6 are shown in figure 7. These loops are 
formed by the lowest two natural frequencies of the beam. The third and 
fourth, fifth, and sixth, and presumably all higher frequencies coalesce 
at larger values of A'. For A' = 0 (no airflow) the frequencies 
behave in the "usual" manner, namely that their values are increased by 
an increase in boundary restraint. Also, an increase in panel stiffness 
(due to increasing f/c) causes an increase in these frequencies. Com- 
parison of figures 7(a) and 7(d) indicates that f/c has a larger 
effect on the clamped beam frequencies than on those of the simply 
supported beam. 

For A ' > 0 the beam vibrational motion induces an aerodynamic 
loading which subjects the beam to a nonconservative force distribution 
(ref. 18). The stability of such a nonconservative elastic system is 
not necessarily increased with an increase in the forces of restraint. 

In reference 23, for example, the destabilizing effect of increasing 
deflectional spring support stiffness at panel boundaries was noted for 
certain panel configurations. In figure 7 this same type of behavior is 
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Figure 7.- Effect of face sheet 
loops; r’ = 2.0, 




(d) f/c= 0.632 

)ending stiffness on frequency 


n = k' = X' 


0 . 



observed for f/c = 0.01 and f/c = 0.20 where an increase in boundary 
restraint from the simply supported (zero moment) to the clamped 
(nonzero moment) edge condition is destabilizing. 

The variation in with r' is shown in figures 8(a) and 8(b) 

for f/c = 0.01 and f/c = 0.632, respectively (n = = X' =0). For 

f/c = 0.01, the clamped edge condition gives a higher flutter boundary 
than the simply supported edge condition when r' is less than about 0.2. 
For 0.2 <r' <10 the simply supported beam is more stable than the clamped 
beam. The flutter boundaries presented in reference 4, for 0^r'_<2 
and T = 0 (f/c =0), are essentially identical to the r’ ^ 2 portion of 
the curves in figure 8(a). In figure 8(b), f/c is equal to 0.632 and 
in this case the clamped boundary is never lower than the simply 
supported boundary, although in the range 0.5 < r'< 0.8, they do 
coincide. At r' = tt^, the curves pass through the solutions obtained 
by Smirnov in reference 7. 

In reference 5, the analysis of reference 4 (t = 0) was used to 
compute flutter boundaries for a wide range of simply supported panel 
configurations. For panels supporting an in-plane tension load in the 
airflow direction the theory predicted that an increase in core shear 
stiffness (decrease in r) could lower This is illustrated in 

figure 9 (taken from ref. 5) which is for n=n=ljky=T=X=0, 
and k^ = -4. (For r = 0, the value kx = -4 represents a tensile load 
with a magnitude equal to the buckling load of the panel at zero dynamic 
pressure.) It is seen that for r = 2 is more than a third 

larger than X^p at r = 0 (infinite shear stiffness) . To determine 
whether the face bending stiffness affects this behavior (decrease in 
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(a) f/c = 0.01 

Figure 8.- Effect of shear flexibility on flutter boundaries of 
sandwich beams; n = k' = x' =0. 







flutter speed with increase in core shear stiffness), the first seven 
frequency loops for the point r = 2, of figure 9, were computed from 
the analysis herein for f/c = 0.01. These results are presented in 
figure 10 where it is seen that the first two bending frequencies (B1 
and B2) coalesce at = 1231. This is the same value for 

obtained with t = 0. Thus, the flutter boundary in figure 9 appears 
to be unaffected by small, but nonzero, values of f/c (or t) . 

B. Effects of Rotary Inertia 

In the preceding section the flutter condition corresponded to 
the coalescence of the lowest two bending mode frequencies as illus- 
trated in figure 10. It was mentioned in the discussion of equa- 
tions (19) and (20) that the inclusion of both rotary inertia and shear 
flexibility in the analysis admits the two additional motions termed the 
thickness- shear and thickness-twist modes. In general, both of these 
motions produce transverse deflections w(x,y,t) perpendicular to the 
panel's undeformed middle surface. In such cases there is the possi- 
bility of frequency coalescence from all three types of motion since a 
transverse deflection from any of the possible modes will influence the 
aerodynamic loading which in turn may change the panel response. 

For a simply supported panel with an isotropic core there is no 
transverse motion associated with the thickness-twist modes (see 
eqs. (81)). In this case, only the bending and thickness-shear modes 
interact with the aerodynamics. The typical variation of the in-vacuo 
(X = 0) values of these bending and thickness-shear frequencies with the 
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rotary inertia parameter, X, is shown in figure 11.^^ The first 14 
bending frequencies, labeled B1 through B14, are nearly independent of 
X whereas the first two thickness-shear frequencies, labeled TSl and 
TS2, are approximately inversely proportional to x. (These bending 
and thickness-shear frequencies were computed from eqs. (20).) 

Except for X, the panel parameters used to compute the frequency 
curves in figure 11 are the same as used in figure 10, where X = 0. 
Thus, the X = 0 values for the bending frequencies of figure 11 corre- 
spond to the intersection of the frequency loops with the A = 0 line 
in figure 10. For X > 0, the frequency loops of figure 10 are altered 
as a result of the thickness -shear modes. This is illustrated in 
figures 12(a) to 12(e). 

The most striking feature of figure 12(a) (X = 0.002) is that the 
fifth bending frequency (B5) has broken away from the sixth bending 
frequency (B6) . The B6 frequency now coalesces with B7 while B5 pre- 
sumably coalesces with some higher frequency mode. Note that for X = 0 
(fig. 10), the first frequency loop determines A^.^ since each succeed- 
ing loop coalesces at successively larger values of A. Of the first 
six frequency loops shown for X = 0.002 (fig. 12(a)) the first loop 
still coalesces at a smaller A than the others. However, each 


^^For panels having identical faces, and cores which are light 
compared to the face sheet weight, the rotary inertia parameter is 
given, approximately, by 
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(a) X = 0.002 

Figure 12.- Effect of rotary inertia on frequency loops of figure 10; 
n = n = 1, r = 2, = -4, = 0, f/c = 0.01; simply supported 
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succeeding loop does not coalesce at successively larger X (note 
B12-B13). This raises the possibility that X^^ may be determined by 
some still higher frequency modes. 

Increasing X to 0.004 results in additional changes in the modes 
which coalesce. As shown in figure 12(b), the B1 mode has now broken 
away from B2. In addition the first thickness-shear natural frequency 
(TSl) is near the B12 frequency, with which it coalesces. Of the seven 
frequency loops shown the fifth loop (formed by BIO and Bll) coalesces 
at a smaller X than the others. An increase in X to 0.006 
(fig. 12(c)) again changes the pattern and the fourth loop (formed by 
B8 and B9) coalesces at the smallest X. 

When X is increased to 0.008 (fig. 12(d)) the B8 and TSl 
frequencies coalesce and become complex at X = 616. Note that this is 
50 percent less than X = 0 (fig. 10). However, when X is 

increased to 709, these frequencies become real again, making the motion 
of the corresponding modes stable. The motion remains stable until a 
second instability region is reached at X = 995. This type of behavior 
(region of stability bounded above and below by unstable regions) was 
also noted in reference 23 for ortho tropic panels mounted on 
deflectional spring supports. 

For X = 0.010 (fig. 12(e)) the in-vacuo values of the Bll and TS2 
frequencies are nearly equal. As X increases to about 320 they become 
equal but remain real (i.e., the corresponding modes vibrate at the same 
frequency but are stable) . The vibrations of these modes remain stable 
until X is increased to 1042 where their frequencies finally coalesce. 
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This frequency crossover behavior, with the frequencies remaining real, 
is similar to results obtained in reference 24 for panels covering an 
enclosed cavity. 

Another effect of the rotary inertia is that it makes the solution 
for Xqx dependent on ky. This differs from the solution of refer- 
ence 4 where the frequency and crossflow loading terms always grouped 
in the characteristic equation as + n^ky, a change in ky merely 
shifting the frequency loops along the axis. This unique grouping 

does not occur for X > 0 (see eqs. (22c), (24), and (27)). 

The sensitivity of the frequency loops to ky for X = 0.010 is 
indicated by comparing figure 12(e) with figure 13. The only difference 
in the parameters selected for these figures is that k^ = 0 in fig- 
ure 12(e) whereas ky = -4 (tension load) in figure 13. One effect of 
the change from ky = 0 to ky = -4 is that the in-vacuo bending fre- 
quencies are increased from to approximately + 4. The in-vacuo 

thickness-shear frequencies, however, are essentially unaffected by this 
change in ky. The frequency loops for both values of ky are quite 
similar, being formed by the same frequency pairs and coalescing at 
corresponding values of A that differ by less than 2 percent. Thus, 
at least for the parameters considered, the crossflow loading term k^ 
has only a minor effect on frequency coalescence. 

C. Effects of an Qrthotropic Core 

One of the effects of an orthotropic core (r^ ^ Ty) is illustrated 
in figure 14, where r^^ ~ 2.0 and ry = 0.5 (all other conditions are the 
same as in fig. 12(d)). It is seen that the thickness-twist frequencies 
(TTl, TT2, and TT3) are no longer independent of A and that 
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Figure 13.- Effect of crossflow loading on frequency loops of 
figure 12(e); n = n = 1, r = 2, kx = -4, ky = 
f/c = 0.01, X = 0.010; simply supported edges. 


-4 





0 100 200 300 400 500 600 700 800 9 



Figure 14.- Illustration of thickness-twist frequency coalescence 
orthotropic core and simply supported edges; n=n=l, r^=2 


ry = 0.5, kx = -4, ky = 0, f/c = 0.01, X = 0.008, p = 0.3. 



they now coalesce with the bending frequencies. For this particular 
choice of parameters it is also seen that the first two thickness-shear 
frequencies (TSl and TS2) do not coalesce with their adjacent bending 
frequencies. Their behavior is now similar to the first bending 
frequency. 

In the remaining examples of panels having orthotropic cores the 
rotary inertia has been neglected (X = 0) and is determined from 

the coalescence of the lowest two bending frequencies. This facilitates 
the calculation of so that the general effects of changes in core 

orthotropy can be determined with much less numerical computation. 

Figure 15 shows the variation in the flutter dynamic pressure 
parameter (X^^ or X^j>) with length-width ratio (n = a/b) for three 
values of the shear flexibility parameter r^ (or r^) and with r^ 

(or r^) held constant. All results are presented in terms of param- 
eters defined by the shortest panel dimension (length or width). Thus, 
for a/b <. 1 the primed parameters (defined by a) must be used whereas 
for a/b >_ 1 the unprimed parameters (defined by b) must be used. The 
solid curves are from the exact solutions derived herein and were com- 
puted with f/c = 0.01 in the equation for x. For the isotropic core 
(rx = r^, r^ = r^) the flutter boundary coincides with the results of 
references 4 and 5 wherein the face sheet bending stiffness was 
neglected (x = f/c = 0) . The flutter boundaries for the orthotropic 
cores (ry. ^ r^) bracket the boundary for the isotropic core except at 
a/b = 0 where the solution is independent of r^. (see eqs. (26), (42), 
and (53)). The dashed curves in figure 15 are from the two-mode 
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Figure 15.- Comparison of flutter boundaries for simply supported 
panels as predicted by two-mode Galerkin solution and exact 
solution: kx=kv=X= 0 , n=l, p=0.3. 



Galerkin solution of reference 9.^^ This approximate solution predicts 
the correct trends for the range of parameters shown but underestimates 
the exact solutions by a considerable amount. 

The frequency loops for the square panel having r^ = 0.4 and 
r^ = 0,2 (peak of solid and dashed boundaries in fig. 15) are shown in 
figure 16. The loop predicted by the two-mode Galerkin solution is 
symmetric about a line midway between the two in-vacuo frequencies 
(X = 0) , whereas the exact loop leans to the right and coalesces at a 
considerably larger value of X. Note that the Galerkin solution 
correctly predicts the in-vacuo frequencies because the two modes used 
in the Galerkin solution are the exact mode shapes for X = 0. 

Figure 17 indicates the variation in X^^ with changes in the 
core shear flexibilities (r^, r^) for a square panel. The solid and 
dashed curves are for clamped and simply supported leading and trailing 
edges, respectively. These curves show the variation in X^^ with r^ 
for r^ = 0.4. The circular symbols locate solutions for isotropic 
cores and indicate changes in X^;.^ with changes in r^ for 
Vy = constant. By noting the changes that occur in X^^ as r^ and r^ 
are varied, from initial equal values to final equal values, one can 
determine the relative importance of r;^: r^,. 

^^In reference 9 this two-mode solution v/as derived for the more 
general case of a panel having both orthotropic bending and orthotropic 
shear stiffnesses. The results shown in figure 15 are for isotropic 
bending stiffness and the corresponding two -mode solution, in terms of 
the notation used herein, is given in appendix D. 
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Figure 16.- Comparison of frequency loops for a simply supported panel 
having an orthotropic core, as predicted by two-mode Galerkin and 


exact solutions; n = n = !> 


0.4, Ty 


0.2, kx — ky — X — 0 , 


. 3 . 



Figure 17.- Variation in with crossflow shear flexibility; n = n = l, 

f/c = 0.01, ii = 0-3; side edges simply supported. 
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For example, consider the simply supported panel. Starting at 
rjf = ry = 0.2, Xgj. equals 325. I£ one follows the path ry = constant 
to r^ = 0.4, Ty = 0.2 to ry = 1.0 (along the r^^ = 0.4 curve), and 
finally ry = constant (i.e., 1.0) to r^ X • 0 ^ ends up at 153, a 

change in of 172 (53 percent). Of this amount, 29 percent is due 

to the change in r^ and 24 percent is due to the change in ry. Thus, 
for a square panel, the two shear flexibility parameters appear to be of 
comparable importance. 

There is apparently no experimental flutter data for sandwich type 
panels available for comparison with the results predicted by the theory 
presented herein. The validity of the theory's representation of 
ortho tropic cores such as honeycomb, however, appears to be established 
by the comparison between experimental and theoretical vibration results 
reported in reference 25. 
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V. CONCLUDING REMARKS 


On the basis of static aerodynamic theory, an exact solution is 
obtained for the supersonic flutter of flat rectangular, biaxially 
stressed sandwich panels. The airflow is assumed to be parallel to 
simply supported side edges. The leading and trailing edges are either 
simply supported or clamped. The mathematical model describing the 
panel motion includes terms which account for rotary inertia, face 
sheet bending deformations, and transverse shear deformations of an 
orthotropic core, such as honeycomb. Numerical results are presented 
in graphical form to illustrate the effects of these physical properties 
on the flutter solution. For the range of parameters considered, the 
following conclusions are noted. 

The face sheet bending stiffness has a negligible effect on flutter 
if the faces are thin compared to the core thickness. In such cases 
flutter boundaries can be adequately predicted by somewhat simpler 
theories that treat the face sheets as membranes. The face sheet bend- 
ing stiffness becomes more important as the face to core thickness ratio 
increases. This is especially true for panels having relatively 
flexible cores. 

The combined influences of rotary inertia and core shear 
flexibility can markedly change the frequency coalescence behavior. 
Failure to account for this combined effect can lead to significant 
overestimates of flutter dynamic pressure values. The rotary inertia 
also causes the flutter solution to be slightly dependent on the 
crossflow in-plane loading. 
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The directional shear stiffness properties of the core are of 
comparable importance for square panels. As the panel width increases, 
the importance of the shear stiffness in the crossflow direction 
decreases. 

A two-mode Galerkin solution gives the correct trends for the 
flutter boundary of a panel having an orthotropic core but significantly 
underestimates the exact solution. Also, an existing exact solution for 
panels having isotropic cores is shown to be correct only for certain 
restricted cases, the solution having been obtained from a differential 
equation and boundary conditions which are incomplete. 

Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, March 1, 1971 
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VI. APPENDIXES 


The following appendixes include explanatory information concerning 
the panel theories of references 14 and 15, the in-vacuo mode shapes, 
and the Galerkin solution of reference 9. 

A. Stiff ness and Inertia Properties 

The sandwich panel theory of reference 14 was developed in terms of 
orthotropic stiffnesses that are defined by the ratios of various 
internal forces and distortions. For example, the bending stiffness 
is defined as the ratio of moment to curvature w when no 

A A. ^ jCX 

Other forces or moments are acting. This approach is applicable to a 
wide variety of panel constructions (e.g., honeycomb core, corrugation 
stiffened) but requires a separate evaluation of the stiffnesses. 

The sandwich panel theory of reference 15 applies to an isotropic 
panel construction for which the form of the face sheet and core defor- 
mations are specified. The strain energy of the structure is expressed 
in terms of the assumed deformations so that differential equations of 
equilibrium and their corresponding boundary conditions can be obtained 
by a variational procedure. The coefficients of the virtual displace- 
ments appearing in the boundary condition equations are expressions 
involving various panel distortions (w for example) and the elastic 
and geometric properties of the core and faces. Physically, these 
coefficients are the boundary forces associated with the corresponding 
virtual displacements. Hence, they can be used to define the panel 
stiffnesses. 

The approach of reference 15 is easily adapted to define the 
stiffnesses for a panel having isotropic face sheets and a core with 
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orthotropic shear moduli (G^ , Gc ) . (The extensional stiffnesses of 

X y 

the core are neglected.) The expressions for these stiffnesses, defined 
in the manner of reference 14, are given in table I. 

The sandwich panel theories of references 14 and 15 do not include 
inertia forces. Inclusion of the transverse inertia loading is straight- 
forward; inclusion of inertia loadings due to longitudinal and rota- 
tional motions requires careful consideration of the panel *s elastic 
deformations, the position of its center of mass, and its elastic axis 
location. 

The panel deformations (as considered in ref. 15) occurring in the 
X, z-plane are illustrated in figure 18. The slope of the panel faces, 
9w/3x, is comprised of a core shear contribution (Yc^j^) an angular 

contribution (otc;^;) • The quantities Yc^ are defined with 

respect to the core orientation given by the line ab. Alternately, 

9w/9x is comprised of the angles Yx and a^, which are defined with 
respect to the line cd drawn between the midpoints of the face sheets. 

The angle y^ is related to the transverse shear stiffness Dq , 
as defined in the manner of reference 14, by the equation Yx = Qx/^Q • 
Thus 



gives the angle in terms of w and Q^/Dq^ which are considered to 

be unknowns in equations (1). From figure 18, is also seen to be 
related to the x-components of the displacements of the face midpoints 
(ui and U 2 ) by the relation 
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TABLE I.- STIFFNESS AND MASS PROPERTIES FOR SANDWICH 
PANELS HAVING HONEYCOMB TYPE CORES 







FACE I 


C 

a 


1 


BEFORE DEFORMATION 



— U 


2 


Figure 18.- Deformed panel geometry. 


tty = 


Ui - U. 


CA2) 


where h is the distance between the midpoints of the two faces. 

The interpretation of Ox and its relation to the elastic axis 
can be determined from figure 19. This figure is a plot of the 
x-components of the displacements shown in figure 18. The midpoint 
displacement of face i is Uj^. It is associated with the axial force 
Nx^ acting on face i. The total axial force is Nx = Nx^ + Nx 2 which 
is used to define an effective displacement 


_ uiEf^fl + U2Ef^f2 

"" " Ef^fl + 

The use of equation (A3) to eliminate either u^ or U 2 from 
equation (A2) gives 

ui =u+dia^ 

U2 = U " 


(A3) 


(A4) 


where 


di = 


Ef2f2 


1 - »Ef^fl -H Ef^fz^ 
Ef^fl 


do - 


yEf^fl + Ef^fa 
di + d 2 = h = c + 


fl + f2 


“x = 


Ui - U2 


w 




Ox 

% } 


(AS) 


The virtual work of the axial forces Ny. is 


6 W = 6 U 2 
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Figure 19.- Elastic axis location and variation of x-component 
of displacement across the panel thickness. 



which, if equations (A4) are used, becomes 


6W = (Nxi + Nx^)6u + (Nx^di - Nx^d2)6a^ 


(A6) 


Thus, the work due to the axial face sheet loads Nx-^ and Nx 2 is equiv- 
alent to the work done by the moment of these forces (about the point 
defined by d^ (or d 2 )) acting through the angle a^^;, plus the work 
done by the total force N;^:^ + Nx 2 acting through the displacement uT. 
Hence, the term d^ (or d 2 ) locates the elastic axis of the panel cross 
section. 

For face i, the virtual work of the inertial forces due to the 
X and y components of the middle surface displacements (u^^ and Vj^, 
respectively) is 


6Wi 


j'i 


pfi^i 


x=o y=0 


3^Uj 


+ 


9^Vj 

8^^ 


6v^Jdx dy 


Inclusion of this term (for both faces) in the analysis of reference 15 
yields the following equation of motion for the rotation of a differ- 
ential element about the elastic axes perpendicular to the xz-plane. 



The left hand side of this equation represents the restoring moment 
about the elastic axis due to the elastic forces. The first term on the 
right hand side is the mass moment of inertia (Iq) of the face sheets 
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about the elastic axis, multiplied by the angular acceleration of the 
face sheets about the elastic axis; hence this term can be interpreted 
as the moment due to rotary inertia. The last term on the right hand 
side arises because the center of mass does not necessarily coincide 
with the elastic axis. The quantity ~ is equal to 

‘^cm(Pf 2 ^l Pf£^2) where d^.jjj is the distance from the elastic axis to 
the center of mass of the two face 'sheets and (pf^fl + is the 

total mass of the two faces. The acceleration 3^u/3t^ is the 
acceleration of the elastic axis in the x-direction. 

It is noted that equation (A7) is a particular case of the general 
vector equation of motion (ref. 26) 

Sp = '^p + pcm X mrp 

where Mp is the moment of applied forces about an arbitrary point p; 
ftp is the time rate of change of angular momentum about p; is the 

vector distance from p to the center of mass of a system of mass m; 
and ?p is the acceleration of the point p. In equation (A7) , point 
p corresponds to the elastic axis. 

The equation of motion for the x-direction is 

-Nx^x ^^y,y " ^ Pf2^2) ^ (u - (A8) 

The term (u - dj.jjjax) is simply the x-coraponent of the center of mass 
displacement. 

It is seen that equations (A7) and (A8) are dynamically coupled 
(i.e., 3^u/9t^ and S^a^/St^ appear in both equations). However, for 
the case of identical face sheets, or faces which satisfy the relation 
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(A9) 


Pfl Pf2 

the elastic axis and center of mass coincide and the equations uncouple. 
In this case equation (A7) reduces to equation (lb) . Hence, the solu- 
tions to equations (la) , (lb) , and (Ic) that are obtained herein are 
strictly applicable to panels having identical face sheets, or face 
sheets satisfying equation (A9) . Also, the stipulation that Nx and 
are constant assumes that the inertia force due to longitudinal 
acceleration of the center of mass is negligible. 

It should also be noted that the inertial forces due to the 
rotational accelerations of the core and faces about their midpoints 
have been neglected. For typical sandwich construction of lightweight 
cores and thin face sheets, these inertia contributions are expected to 
be negligible compared to the inertia of the face sheets about the 
elastic axis. 

B. Comments on So lution of Reference 27 
It was noted on page 26 that the thickness-shear and thickness-twist 
motions are due to the combined influence of core shear deformations 
and rotary inertia moments. In reference 27 the flutter analysis of 
reference 5 (isotropic core, face bending stiffness and rotary inertia 
neglected) was extended to include the effects of rotary inertia. 
Unfortunately, the analysis of reference 27 is based on differential 
equations that do not properly account for the rotary inertia moments. 

The moments due to rotary inertia are proportional to the 
acceleration of the rotation of the panel cross sections, for example, 

= [w,^ - (Qx/Dq^)] (see appendix A). In reference 27, however. 
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the rotary inertia moments are taken proportional to the acceleration of 
the sZoipee of the panel cross sections, for example, The 

rotation and the slope of a cross section are not equal unless the 
core shear deformation vanishes (Dq = «>) ; hence the analysis in 
reference 27 leads to erroneous results since it is based on an 
incorrect description of the rotary inertia moments. For example, the 
characteristic equation of reference 27 is linear, rather than cubic, 
in 0 )^ and does not predict the thickness-shear or the thickness-twist 
frequencies. Also, the "bending” frequencies are correct only if the 
rotary inertia or the shear flexibility is neglected. Consequently, the 
results presented in reference 27 are of questionable value. 


C. In-vacuo Bending and Thickness-Shear Mode Shapes 
for a Simply Supported Sandwich Beam 
The two differential equations governing the beam behavior (n = 0) 
are obtained from equation (82) . 


Lii Li2 

J-21 L 22 

The differential operators Lij are given by equations (83) 
that all terms involving derivatives with respect to y are 
For simply supported edges at x = 0 and x = a the boundary 
are 



(Cl) 


except 

omitted. 

conditions 
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w = 0 




w =0 
,xx 


(C2) 


Expressions for the transverse deflection w and the shear angle 
that satisfy the boundary conditions are 


, . V . . im:x iu)t 

w(x,t) = Am sin e 


YjjCx,t) = 


Qx(x,t) Bm 

mTTX icot 


fC3) 




^ COS — — e 

% 


The above expressions for w and Yx also satisfy the differential 
equations, provided that 


-Tin^+m^k'+u' ^ 


-m 


m 


(-m^+x'w' 1-r ' (-m^+x'oj'^) 



/ 

] 





> 

> = 

.0, 


.tt^Ds 



CC4) 


For the nontrivial case, equation (C4) is satisfied by equating the 

determinant of the square matrix to zero. The resulting equation can be 
_ 2 

solved for to' and yields the following frequency equation: 


-I 2 

“b 

TS 


2 [m‘++ (Tm‘*-m2ky) (1+r^m^) ] 


[1+ (ri+X ' )m2+r^x ' (^‘^-m^kj;.) ] ± J [l+Cr^^-X ' )ni2-rxX ' (Tm''-m2k^) ] ^+4X 


'm2 


CC5) 


As in the case of the panel (n > 0, eq. (20a)) the smaller solution for 
(o' (+/~ ) gives the frequencies of the bending modes (tOg ) , and the 

larger solution (-/~ ) gives the frequencies of the thickness-shear 
modes (m.j,g) • 
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The shapes o£ the bending and thickness -shear modes are obtained 
from the relation between Am and This relation, determined from 
equation (C4) , is 



1 

m 


1 + r^^ m2 


V ' — ' ^ 

X o)g 
TS/ 


4K 


X’oT' 


2 

B 

TS; 


or 



m 


+ m^kv + o) 


I 

•XI 


I 2 

X + 

TS 


(C6) 


When equations (C6) are used in conjunction with equations (C3) , the 
following relation between the slope of the beam (w and the shear 
angle (y^^) is obtained 


r' (m2 - x'oo 


_t 2 

B 

TS 


w 


,X 


1 + r' (m2 - x'tOr^ 


“B 

TS, 


(C7) 


or 


X ^X (—'2 2t ' 
= — (“b ® ^X 

TS 


w 


xm^ 


,x m'’ 


It is recalled from appendix A that 


W = tty + Y„ 

,x ^ 'X 


(C8) 
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where is the beam rotation about the elastic axis. The transverse 

deflection is thus 

W = /(«x + 

and is comprised of a contribution due to rotation and a contribution 
due to shear. The proportion of each contribution is determined as 
follows: 



where equations (C7) and (C8) have been used. Thus, the relative 
contribution to w (or w produced by the rotation is given by the 
first term in the bracket of equation (C9) ; the second term gives the 
relative contribution due to shear. Note that the shear contribution 
vanishes for r^ = 0 (infinite shear stiffness) . Additional information 
about the shear and rotational contributions is revealed by writing 
equation (C5) in the following form. 


TS 


2 [m^+ 


[l+(r^+X')ni2+r^X» (im^-m^k )] (l ± I 1 - — 

' [U(r;+X')m2-Hr;x’(Tm^-m2k;)]M 


(CIO) 
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For many practical parameter values^ ^ the magnitude of the underlined 
fraction in equation (CIO) has a value much less than 1.^^ In such 
cases equation (CIO) can be approximated by 


“b 


2 , 
TS 


2 (1+rxtn^)] 


[l+(rx+X (im^-m^kx)] 1 ± < 1- 


2rxX* [m^+(Tm^-m^k^) (1+rxm^)] 
[1+ (r^+X • )m2+r^X ' (rm^'-m^kx) ] ^ 


(Cii) 


which in turn leads to 


2 1 + (r^^ + X')m^ 

(ol,^ ^ 

TS 


r^X- 


+ Tm^ - m^k^ 


(C12) 


and 




m IT + i—t) - 

1 + r^m^y X 

A 

1 + a 


i2k: 


(C13) 


For identical face sheets and f/c £ 0.632, 

^ Y 

^ 0.05 


1 

" = 3 


.1 + 


Positive values of the in-plane load parameter are limited by the 
buckling value, k^ = t + [1/(1 + r^^)] (for no airflow). For lightweight 
cores and identical face sheets the ratio X'/^x is approximately equal 
to 




Numerically, this ratio often has a value considerably less than 1. 
^^One notable exception is when r^tm^ >> X'xm^ ^ 1. 
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where 


a 


m^x* 



Th I f 

-f r^-m 

1 + 



In many instances o is also much less than 1. For example, i£ 

= X = 0 then a = m^X'/(l + < ^'/^x which is normally small 

compared to 1. 

Substitution of the above approximate expressions for the 
thickness -shear and bending frequencies into the second of equations CC7) 
yields 



1 + 


1 + r^m^ ^ 

m2x* 


> 



1 


9 * 

2 * 

+ m^r^ 


/ 


(C14) 


where a has been neglected with respect to unity. Since w^jj = ax + Yjf, 
and since no term on the right hand side of equations (C14) can be less 
than zero, the approximate bounds on and are 

Yx 


0 :< 


w 


< 1 




0 < 


a-si 


w 


< 1 




> for the bending modes 


(CIS) 
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for the thickness -shear modes 


(C16) 



The mode shapes resulting from these restrictions on Yx 

illustrated in figures 3(c) and 3(d1 for the bending and thickness -shear 

modes, respectively. 

2 

Niamerical values for w' and Y^/w ^ predicted by the 
approximate expressions (eqs. (C12) , (C13) , and (C14)) and by the exact 
equations ((C5) and (C7)) are compared in table II for three sets of 
parameter values. 
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TABLE II.- COMPARISON OF IN-VACUO BENDING AND THICKNESS-SHEAR 


FREQUENCIES AND MODE SHAPES AS PREDICTED BY EXACT AND 


APPROXIMATE SOLUTIONS FOR A SIMPLY SUPPORTED BEAM 






1 

2 

3 





0.0 

-1.0 

0.0 



f/C 

0.04 

0,04 

0.0 



X' 


0.01 

0.01 

0.002 



f 

’^x 


1.0 

0.1 

0.2 



m 


2 

1 

2 




Exact 

3.203 

1.894 


_,2 

0) 

Bending 


Approximate 

3.208 

1.910 


Thickness- shear 


Exact 

500.8 

1109 






Approximate 

504.0 

1111 





Exact 

0,799 

0.0894 

0.443 

Yx 

Bending 


Approximate 

0.800 

0.0909 

0.444 


Thickness -shear 


Exact 

125.2 

110.8 

225.6 




Approximate 

126.0 

111.0 

226.0 
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D. Two-Mode Galerkin Solution of Refereiice 9 
In reference 9 a two-mode Galerkin solution for flutter was 
obtained for a simply supported panel having: 

1. Two transverse shear stiffnesses (Dq^, Dq ), 

2. Two bending stiffnesses Dy) . 

3. Two Poisson ratios Py) • 

4. A twisting stiffness Dxy 

For isotropic face sheets and an orthotropic core with negligible 
bending stiffnesses, one can set (ref. 14) 

= t'y = ^ 

Dx = Dy = Dg (1 - p2) 

^xy “ ^s 

By use of the above relations, the two-mode solution (in terms of the 
notation used herein) is given by 


A = |i^ - j2|/-p(i)p(j) 


(Dl) 


where i and j are integers whose sum is an odd integer and 


p(m) = 


[-c2Cm) - c3(m) . c4(m)] 


(a) cl(m) 


2 


cl(m) = 1 + 


(^)te) *(¥■) Vy[(") * [(") 


cl(m)kj^ - 2n2 


c2(m) 

c3 (m) 
c4 (m) 


cl(m)[o)2 + n^ky] - n^^ 
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The solution for (at BX/dui^ = 0) is given by 


^cr = fij - J^|[E(i) - E(j) + n'+FCiJ)] 


(D2) 


where 


E(m) = 


c2(m) + c3Cm) 


cl (m) 


and 


F(i,j) 


^ cl(j) - cl(i) 
cl(i)cl(j) 


Equations (Dl) and (D2) provide an estimate for the flutter 
behavior of a simply supported sandwich panel having an orthotropic 
core and isotropic face sheets. Note that this solution does not 
account for rotary inertia, nor for the face sheet bending stiffness 
contribution to lateral equilibrium (i.e., X = t = 0) . 

For X = 0, equation (Dl) can be solved for and, except for 

the absence of t, yields the same in-vacuo frequency equation as given 
by equation (18a). This observation allows equation (Dl) to be written 
as 



where e = - [(mj^ + ijo^)/2] is shown in figure 20. The quantity oT? 

is the in-vacuo frequency of the ith mode (i.e., m = i in eq. (18a)). 
It is seen that X(aT^) is symmetric about the line (oT? + 00^/2. 
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